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Abstract 

We construct a family of random Jordan curves in the plane by 
welding together two disks on their boundaries using a random homeo- 
morphism. This homeomorphism arises from a random measure whose 
density, in a generalized sense, is the exponentiated Gaussian Free Field 
at criticality. We also introduce a representation of the Gaussian Free 
Field in terms of vaguelets, which may be of separate interest. 

The result extends a theorem of Astala, Jones, Kupiainen and Sales- 
man f |AJKS09j ) to criticality. 

1 Introduction 

We extend a theorem of Astala, Jones, Kupiainen and Salesman QAJKS09J) 
to criticality. The authors of [AJKS09] constructed a family of random Jor- 
dan curves in the plane by solving the conformal welding problem with a 
random homeomorphism. The random homeomorphism arises from a ran- 
dom measure which is constructed by a limiting process via the exponenti- 
ated Gaussian Free Field (its restriction to the unit circle). The algorithm 
depends on one parameter, inverse temperature, and provides a random 
Jordan curve for each inverse temperature less than a certain critical value. 

In this work we extend the construction to criticality. While Astala, 
Jones, Kupiainen and Saksman used a white noise representation for the 
Gaussian Free Field, we use a vaguelet representation. 

Over the last decades there has been an interest in conformally invariant 
fractals which could arise as scaling limits of discrete random processes in 
the plane. We provide an instance of such fractals in the form of Jordan 
(simple, closed, locally connected) curves. 

One of the most important example of conformally invariant fractals 
is Schramm-Loewner Evolution (introduced by Schramm, [SchOO]). One 



version describes a random curve evolving in the disk from a point on the 
boundary to an interior point. The second version describes a curve evolving 
in the upper half plane from one point on the boundary to another point on 
the boundary (typically infinity). The construction depends on a parameter 

[0,8]. It has been shown that the scaling limits of some discrete processes 
are SLE K curves. For example, Lawler, Schramm and Werner ([LSW04J) 
proved that the loop erased random walk converges to SLE2. Another 
example is the percolation exploration process which converges to SLEq on 
a triangular lattice as proven by Smirnov ( |Sm01| ). For more information on 
SLE K and an introduction to percolation we refer the reader to G. Lawler's 
book ([LU5]) and W. Werner's notes ( |W09] ). 

Unlike SLE K , the random curves constructed in the present work (and 
in [AJKS09J) do not evolve in the plane and are closed. Our curves come 
in two varieties. The first is the result of welding a deterministic disk to 
a second one by a random homeomorphism which arises by exponentiating 
the Gaussian Free Field. The second variety is the result of welding two 
disks on their boundaries by a random homemorphism which arises from 
two independent Gaussian Free Fields. 

The subcritical case (inverse temperature less than the critical value) has 
also been studied by Sheffield ( [S10J ) . Using different methods he welded 
together two disks using two independent Gaussian Free Fields to obtain 
SLE K (where k equals twice inverse temperature). Starting with a GFF on 
a disk, Sheffield defines a random area measure on the disk. The conformal 
equivalence class of the disk and random area measure is a model for ran- 
dom surfaces known as Liouville quantum gravity. The surfaces are called 
quantum surfaces. Sheffield welded together two quantum surfaces of nor- 
malized quantum area by matching quantum length on the boundaries. The 
two GFF version of Astala, Jones, Kupiaien, Saksman welds together two 
quantum surfaces by matching normalized quantum length on their bound- 
aries. It is not clear whether the two GFF version of [AJKS09] are SLE K 
curves. Binder and Smirnov claim that the two sided version of the subcrit- 
ical Jordan curves look like SLE K when the local dimension of the curves is 
considered ( [BS12] ). In addition, they proved that the one sided version is 
not SLE K (reported by Sheffield in [S10J). 

To construct the closed random curves we follow the broad framework 
developed in [AJKS09]. We define a random homeomorphism by expo- 
nentiating (and normalizing) the Gaussian Free Field. We then solve the 
conformal welding problem for this homeomorphism. While the approach is 
the same as in [AJKS09J, the details are subtantially different. 

The Gaussian Free Field is a random distribution of great importance 
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in statistical physics and has been studied extensively in both the physics 
and the mathematics literature. For a mathematical introduction see [S07J 
as well as the introductions of |SS0 9] and jSSlOj . In [AJKS09], the authors 
used a white noise representation for the GFF. We use a vaguelet represen- 
tation which was suggested to us by Peter W. Jones. Vaguelets are functions 
very similar and related to wavelets and appear, for example, in the work of 
Donoho ([D95]) and Meyer and Coifman QMC97J). The vaguelet represen- 
tation allows us to work with the Gaussian Free Field scale by scale. 

The procedure of exponentiating the Gaussian Free Field to get a random 
measure and then a homeomorphism appears also in the work of Sheffield 
QS10J) and Duplantier and Sheffield(e.g. jDSllj ). The measure constructed 
is a type of multifractal multiplicative cascade as introduced by Mandelbrot 
( |M74j ) and studied, among others, by Kahane ( |K85j ). Kahane and Peyriere 
( [KP7H| ). Bacry and Muzy(jBMQ3]) and Robert and Vargas ( |SV08] ). Mul- 
tiplicative cascades appear also in mathematical finance as an important 
part of the Multifractal Model for Asset Returns introduced by Mandel- 
brot, Fisher and Calvet ( |MFC97 J). 

Most of the effort in this work is spent on solving the conformal welding 
problem. The classical result on this topic is that quasi-symmetric homeo- 
morphisms are welding maps (see e.g. jLV73] ) and thus give rise to conformal 
weldings. Lehto ([L70]) and David ( [ D88 ] ) relaxed this assumption. Oikawa 
( |U61| ) and Vainio ( |V85| . [V89J and |V95| ) provided examples of home- 
omorphisms which are not welding maps. Hamilton relaxed the concept 
of welding map (introducing the generalized welding) in [H91J and Bishop 
proved that every homeomorphism is almost a conformal welding in the 
sense that one can modify it on a small set and make it a welding map (see 
|B06| ). We refer the reader to [H91J and [B06J for more information on the 
welding problem and related topics. In our setting the random homeomor- 
phism fails to be quasi-symmetric and we have to employ tools developed 
by Lehto in [L70J to get a welding. 

We also prove that the random curves we construct are unique up to 
Moebius transformations. The uniqueness is equivalent to the uniqueness 
of the solutions of the Beltrami equation (see section [3] for more details). 
While for quasi-symmetric homeomorphisms this follows from the measur- 
able Riemann mapping theorem, in our case it is a consequence of a deep 
theorem by Jones and Smirnov ([JSOOJ) and its extension by Nienminen 
and Koskela ([KN05J) on conformal removability. There are cases of "wild" 
homeomorphisms where the welding is far from unique (see [B94J and [B06]). 

The present work proceeds as follows: In section [2] we state the main 
theorems. In section [3] we introduce the conformal welding problem and the 
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approach we use to solve it. We also state the main probabilistic estimate 
we need to prove the theorems and solve the conformal welding problem. 
Section [4] presents the construction of the vaguelets and proves some of their 
properties. Section [5] introduces the Gaussian Free Field and its representa- 
tion in terms of Fourier series, vaguelets and white noise. Section [6] outlines 
the construction of the random measure, its properties and presents also the 
proofs of these properties. The last part gives a proof of the main probabilis- 
tic estimate by describing the decoupling and estimating the distributional 
properties of the random variables involved (section [7]). It then describes the 
construction of the random tree and its survival properties (section [8]), gives 
the proof of the main modulus estimate we need (section |9| and completes 
the probabilistic estimates. 

Acknowledgments The present work was done during my graduate 
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Jones, for proposing this problem and for his support during my PhD. I am 
also grateful to Ilia Binder for several discussions on the topic. 



2 Results 

Following [AJKSQ9] we can define the restriction of the (random distribu- 
tion) Gaussian Free Field on the circle by: 

A n cos(2irnQ) + B n sin(27rn#) , . 

— ^ , In 



Jn 

71=1 



where A n ,B n ~ -/V(0, 1) are independent. In what follows we will use the 
alternative representation 

i 

where Ai ~ iV(0, 1) are independent and {tpi} are periodized half-integrals 
of wavelets. They are called vaguelets and appear for example in Donoho 
([D95]) and Meyer and Coifman QMC97J). Vaguelets satisfy essentially the 
same properties as wavelets. 

We can now consider a sequence t n and define the following random 
measure 

i,,= ^e^- k{ '-* m -^' {,m m (2) 

where ai ~ -/V(0, t log i ). Kahane proved in |K85j that if t n = t < t c = 2 

this limit exists, is a non-zero, finite, non-atomic singular measure almost 
surely. 
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We extend the result to certain increasing sequences t n — > t c = 2. More 
precisely, given 7 < 1 we take t n = t c — kr 1 for all n G [nfc,nfc +1 ), where 

n k ~ {k + l)Te c '^ c+1 ) 37(fe+1) / e , C is a large constant and e > small. 
The next step is to define the random homeomorphism h : T — > T: 

fe(0) := i/([O,0))/i/([O,2tt)) for G [0,2tt). (3) 

Our goal is to solve the conformal welding problem for this homeomorphism: 
we seek two Riemann mappings /+ : D — > and /_ : Boo — > $7_ onto the 
complement of a Jordan curve T such that h = (f+)~ 1 /-■ 

Theorem 1. Almost surely, formula^ defines a continuous circle homeo- 
morphism, such that the welding problem has a solution T. The curve T is 
a Jordan curve and bounds a domain £1 = f+(B) with Riemann mapping f + 
having the modulus of continuity better than 5(t) = e~^ log * - 1 . For a given 
realization uj, the solution is unique up to Moebius transformations. 

As mentioned above, this theorem is an extension of the result of Astala, 
Jones, Kupiainen, Saksman QAJKS09J. In their paper, the variances t n 
were all equal to a value strictly less than the critical one. In addition, the 
Riemann mapping /+ was Hoelder continuous. By contrast, the variances 
here tend to the critical value, and the Riemann mapping satisfies a weaker 
modulus of continuity. 

The construction presented here can be used to construct weldings also 
for sequences tk which converge to t c = 2 faster than in the theorem. How- 
ever, the modulus of continuity that we obtain in those cases will be worse 

than e V^ log * and we cannot prove the uniqueness of the welding up to 
Moebius transformations. Uniqueness, in our setting, is a consequence of 
the removability results of Jones and Smirnov QJSOOJ) and Koskela and 
Nieminen ( KM)')! ) for conformal mappings. These results provide a suffi- 
cient modulus of continuity for conformal removability, which, in turn, im- 
plies uniqueness of the welding. It is not known what the optimal modulus 
of continuity is that ensures removability (see the papers [JSOO], [KN05j, as 
well as the paper by Jones and Makarov on harmonic measure [JM95J). 

Consider now two independent Gaussian Free Fields and construct two 
independent random homeomorphisms hi and hi- The proof of the following 
theorem is the same as of the first. 

Theorem 2. Almost surely, formula^ for two independent GFFs defines 
continuous circle homeomorphisms, such that the welding problem for home- 
omorphism h\ o h^ 1 has a solution T. The curve V is a Jordan curve and 
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bounds a domain f2 = /+(D) with Riemann mapping f + having the modulus 
of continuity better than 5{t) = e~^ log . For a given realization oj, the 
solution is unique up to Moebius transformations. 

The second theorem is related to a result of Sheffield ( [SlOj ) . He proved 
that welding two disks using subcritical GFFs yields SLE K with k = 2t. 
However, there are two differences. The first difference is that he welds two 
disks after having normalized their random area (defined in a similar way 
to the measure v above). In this construction, as well as in AJKS09 , one 
first normalizes the random length of the boundary. The second difference 
is that Sheffield deals with the subcritical case and we deal with the criti- 
cal one. One can ask whether our construction produces curves which are 
SLE^. The answer is probably no. If the Astala, Jones, Kupiainen, Saksman 
construction yields SLE K , our curves are perturbations of different SLE K 
curves at different scales. However, as Peter Jones has pointed out in a 
private communication, these do not converge to SLE4. 

Astala, Jones, Kupiainen and Saksman prove their theorem by first prov- 
ing estimates on the random measure and then solving a degenerate Beltrami 
equation using a criterion of Lehto (see section [3] for references and details) . 
In order to apply this criterion they have to decouple the distortion and 
prove that the different scales behave roughly independently. 

In the present work we combine the properties of the vaguelets with a 
martingale square function argument to prove the existence of the random 
measure v. We also describe several of its properties. We then do a decou- 
pling, a stopping time argument and a modulus estimate to ensure that with 
high proability the distortion in the Beltrami equation does not diverge too 
fast. Finally, one has to do a careful decoupling of the Gaussian Free Field 
via vaguelets to show that the scales are practically independent. 

3 The welding problem 

Our goal is to contruct a random Jordan curve in the plane. We accomplish 
this by solving the conformal welding problem with a random homeomor- 
phism. Let T, D and Doo be the unit circle, open unit disk and the comple- 
ment of the closed unit disk respectively. The conformal welding problem is 
as follows. 

Let 4> be a homeomorphism on T. We seek two Riemann mappings 
/+ : D — > 0+ and /_ : D ro — > f2_ onto the complement of a Jordan curve T 
such that 4> = /-■ 



6 



Our strategy is to find these mappings by solving the Beltrami equation. 
Assume that <j> = /| T where / G W£r(D,D) n C(B) is a homeomorphism 
and a solution of the Beltrami differential equation 



df 

dz 



H{z) 



df 

dz 



z £ n = d. 



p is called Beltrami coefficient. If the following Beltrami equation 



OF 

dz 



l (z)p(z) 



OF 

dz 



zeC. 



has a unique (normalized) solution F, we can take /_ = -FIdoo- The unique- 
ness implies F can be factored as F = /+ o / for a conformal mapping 
i^(D). Then we will have <f> = (Z^)^ 1 o /_. The road we take is 
now plain: starting with the homeomorphism (ft, we find a mapping / which 
is an extension of (ft to the disk. Then we solve the Beltrami equation on 
the entire plane. Two questions arise: how do we extend the function (ft to 
the disk? and can we solve the Beltrami equation for that extension? 

It is part of classical complex function theory that the Beltrami equa- 
tion admits unique (normalized) solutions whenever ||/x||oo < 1 (m is called 
uniformly elliptic). An extensive reference on the topic is [AIM09]. 

The homeomorphic solutions (in W 1,2 (C)) of the uniformly elliptic Bel- 
trami equation are called quasiconformal mappings and have many inter- 
esting properties. While a conformal mapping maps infinitesimal disks to 
infinitesimal disks, quasiconformal mappings map infinitesimal disks to in- 
finitesimal ellipses. The Beltrami coefficient p is also called "ellipse field" 
and describes this infinitesimal correspondence. 

An important quantity associated with a quasiconformal mapping/Beltrami 
coefficient is denoted by K = j^j^j and is called the distortion of the map- 
ping. If ll/ulloo < 1, the distortion is bounded. 

An equivalent definition of quasiconformality is in terms of moduli of 
annuli. The modulus of an annulus can be defined in several equivalent 
ways. We first introduce the modulus of a family of curves. Given a family 
of locally rectifiable curves V in domain C C the (conformal) modulus is 
the quantity 



where the infimum is over all metrics p : f2 — > [0, oo] such that J pds > 1 for 
all curves 7 C V. Such metrics are called admissible. 



mot 



>d(T) = mod n (T) = inf{ / p 2 dA(x,y)} 



(4) 
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The modulus is a conformal invariant: mod(T) = mod(F(T)) for any 
conformal F. For quasiconformal mappings with distortion K we have 
mod(F(T))/K < mod(T) < Kmod(T) for any family of curves T. In fact, the 
previous definition is equivalent to the almost invariance of the conformal 
modulus (for more details see the book [LV73| ). 

The modulus of an annulus is given by mod(A(r, R)) = ^log(^) and 
it is the same as the modulus of the family of curves which separates the 
two pieces of the complement of A. A third equivalent definition of quasi- 
conformality is in terms of moduli of topological annuli: a mappings F is 
quasiconformal if and only if the modulus of topological annuli is distorted 
by at most a multiplicative factor K. 

The mappings <j) which can be written as /|t with fi uniformly elliptic 
are called quasisymmetric and satisfy K = sup s t < oo. 

In our setting (and the one of Astala, Jones, Kupiaien, Saksman), the 
mapping 4> is not quasisymmetric. While it can be written as the restriction 
of a mapping /, the corresponding fj, is not uniformly elliptic. The mapping 
/ is called degenerate quasiconformal because the distortion, while finite 
almost everywhere, is unbounded. 

O. Lehto ([L70]) proved a very general criterion which ensures the exis- 
tence of solutions to the Beltrami equation in the degenerate case. 

Define: 

1 dp 

J™ K{z + pe id )deJ 

Theorem 3 ( [AIM09] . p. 584). Suppose \i is measurable, compactly sup- 
ported and satisfies \p(z)\ < 1 almost everywhere on C. If the distortion 
function K = jzj^\ * s locally integreable and for some Rq, the Lehto integral 
satisfies: 

L K (z,0,R ) = oo,Vz G C (6) 
then the Beltrami equation 

df , ,df 

- = p(z)-,zec. 

has a homeomorphic solution in W,'. 

We will be using a version of this theorem to prove the main result. 
Lehto's result says that if around every point one can find an infinite number 
of annuli whose conformal modulus is not distorted much by a mapping with 
Beltrami coefficient p, then the Beltrami equation has a solution. We will 



L K (z,r,R) = I — g (5) 
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prove that we can indeed find an infinite number of "good" annuli (i.e. 
annuli which are not distorted too much) around every point. 

Astala, Jones, Kupiainen, Saksman worked directly with the Lehto inte- 
gral above. We will work with moduli of annuli. 

To solve the welding problem (i.e. prove the existence) it suffices to have 
the local uniqueness of the solution of the Beltrami equation. This we have 
because the distortion degenerates as we move closer to T = dD and it is 
bounded inside the unit disk. 

The uniqueness of the welding up to Moebius transformations is equiv- 
alent to the global uniqueness of the solution of the Beltrami equation. 
The information we have is not enough to draw this conclusion immedi- 
ately. We have to apply a deep theorem of Jones, Smirnov( [JS00j ) /Koskela, 
Nieminen([KN05j) on conformally removable curves. 

A curve V is called conformally removable if every global homeomorphism 
which is conformal off T is automatically conformal on C. If we knew the 
curves we contruct were conformally removable then we would automatically 
get the uniqueness of the welding, as any two weldings give rise to a global 
homeomorphism which is conformal off the curve T. 

Jones, Smirnov([JS00j) /Koskela, Nieminen(|KN05]) gave sufficient con- 
ditions on the modulus of continuity on the Riemann mapping /+ that ensure 
conformal removability. As long as the modulus of continuity of /+ is better 

than t — >■ e c V log * (for a large constant C), then the curve is removable. 

We stated in section p that we can solve the welding problem with home- 
omorphism cj) = hi o h>2 • The reader is encouraged to think about this as 
follows (in this way the proof will be exactly the same as for one GFF). For 
each hi consider extension /j. For h\ the extension is to the unit disk D. 
For /12 the extension is to the complement of the unit disk Dqo . Each exten- 
sion fi has a Beltrami coefficient Hi. We seek a solution F to the Beltrami 
equation with coefficient: 



By the (local) uniqueness of the solutions to the Beltrami equation there 
is a conformal map /+ : D — > F(D) such that f = /+ ° /i on D and a 
conformal map /_ : Dqo — > i^D^) such that F = /_ o f 2 on . Then we 
have /+ o h\ = /_ o /i 2 on dD. 
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3.1 The welding problem in our setting 

In section [6] we construct a random measure and we prove, among other 
properties, that it is almost surely not zero on any interval. This allows 
us to define a random homeomorphism <j> on T. Our goal is to solve the 
conformal welding problem for this homeomorphism. 

As espressed in the introduction the broad framework is the one from 
|A JKS09] . Let h(x) = ^jjjyxjy, x £ [0, 1) where v is the measure we construct 
in [6j Extend h periodically to R by setting h(x + 1) = h(x) + 1. Extend 
h to the upper half plane by setting (following Ahlfors-Beurling; see e.g. 
|AIM09j ) 



F(x + iy) = - [ (h(x+ty) + h(x-ty))dt+i f 
2 Jo Jo 



(h(x+ty)-h(x-ty))dt (7) 



for < y < 1. This function equals h on the real axis and it is a continuously 
differentiable homeomorphism. For 1 < y < 2 define F(z) = z + (2 — y)co, 
where Co = Jq 1 h(t)dt — 1/2. For y > 2 define F(z) = z. We also have F(z + 
k) = F{z) + k. On the unit circle we define the random homeomorphism: 

(p(e 2nix ) = e 27Tih ( x \ 

and the mapping: 

<£(» = exp(27riF(logz/27ri)),z £ D 
is the extension of <f> to the disk. The distortions of F and ^ are related by 
K(z, *) = F), z = e 2 ^, to G M+. 
In addition define 

^) = ^-^D; /i(z)=0,z^D. 

Our goal is to solve the Beltrami equation with this Beltrami coefficient. 

We give now an upper bound for the distortion of \i that we will use. 
We also introduce necessary notation. 

Let T> n be the collection of dyadic intervals of size 2~ n . For a dyadic 
interval /, let be the union of / and it's to neighbors of the same size. Set 
d = {(x,y)\x G /, 2-"- 1 < y < 2~ n }. Following |A.TKS09j let J = {Ji, J 2 } 
and set 
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In addition, define 

J(I) = {J = (Ji, J 2 ) : Ji E V n+5 , Ji C (9) 

and 

k u (i)= E <^ J ) ( 10 ) 

Jcj(/) 

The distortion of ^ is the same as the distortion of F (the points are 
mapped appropriately). In the upper half of the square with base /(denoted 
by Ci) the distortion of F is bounded by CqK v (I), for a universal constant 
Co- As a consequence, studying the distortion of p is really about studying 
the doubling properties of the random measure v in In the rest of this 
paper we will only use K v . 

3.2 Main probabilistic estimate 

We want to prove that almost surely we can find infinitely many annuli 
around each point on the unit circle which are not distorted much by a 
mapping with Beltrami coefficient p. We will need the following theorem: 

Theorem 4. There are sequences p n , p n , N n ,b n ,c n such that: 

N n 

P(^2 Mod(G(A(z, p n pl 2PM)) < c n N n ) < p nP £ +K)N - (11) 

i=l 

for any z E T and any mapping G with Beltrami coefficient p. 

We apply this theorem to a net of points on the unit circle and use 
the Borel-Cantelli theorem to get the desired statement that almost surely 
around every point on the unit circle there are infinitely many annuli which 
are not distorted by much. 

This result replaces the Lehto estimate (theorem 4.1) from Astala, Jones, 
Kupiainen, SaksmanQAJKS09]). While their estimate covered scales one to 
p n , this estimate deals with the scales in chunks. 

The theorem is a statement about distortion. Ideally the distortion in 
one scale would be independent of the distortion in another scale. How- 
ever, this is not the case, each scale being correlated with every other scale. 
Fortunately, the correlations decay exponentially. 

The setting here is more complicated than the one in |AJKS09"] . They 
used a representation of the Gaussian Free Field in terms of white noise W 
which had a simpler correlation structure. 
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3.3 Solution of the random welding problem 



rl,l 



Theorem 5. Almost surely there exists a random homeomorphic W,' — 
solution f : C — > C to the Beltrami equation dgf = pd z f , which satisfies 
f(z) = z + o(l) as z — > 00 and whose restriction to T has the modulus of 
continuity uj(t) < e~P° s *^ 

Proof. The proof is essentially the same as in [AJKS09] . In this proof we 
use estimate (11) 

We start by considering for each n an [p^Pn ^ 1+bn / 2 ^ Nn j =: Tn _ ne ^ Q f 
points on [0,1] and denote Cn,fc = exp(27ri/c/r n ) for k G {l,...,r n }. Set 
also G n = {Cn,l; • ■ ■ )Cn,r„}- Any other point on T is at distance at most 
~ Pnpn +bn Nn horn G n . Define the event 

A n , k := {Mod(F{A(Q n ^pnP^ +K/2)N \Pn))) < c n N n ] 

Now set A n = C\kA n ^. Since 

E^n) < EE^) * E^^ B * E^ W2 * E ^72 < ~- 

n=l n fc=l n n n 

Borel-Cantelli tells us that almost every realization u is in the complement 

of U n >n (w)An- 

Consider the approximations p\ = A^p to p. For each I denote by fi 
the normalized (random) solution of the Beltrami equation with coefficient 
pi and such that fi(z) = z + o(l) as 2 — >• 00. In other words /; is a quasi- 
conformal homeomorphism of C. This solution is obtained by means of the 
measurable Riemann mapping theorem (sec [A06J ) 

We want to prove that almost surely the family {/;} is equicontinuous. 
Outside D all these mappings are conformal and equicontinuity follows from 
Koebe's theorem. Equicontinuity inside D follows from the fact that at any 
point inside the disk the distortion is determined by the measure v on finitely 
many intervals and hence it is bounded. 

To prove equicontinuity on T we consider the functions Fi(z) = fi(e 2mz ) 

For P-a.e. oj we have 

n 

Mod(F l {A(C n , k ,Pnp£ +Kl2)N \l)))> E d^Vl 

i>i (£j) 

Lemma 2.3 in |AJKS09| gives: 
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for any F quasiconformal. 

Fix one realization oj. Putting together the last two inequalities (Fi is 
quasiconformal for any I). 

dmmiF^BiC^pnp^/ 2 ^))) < 16diam(FKi?(Cn,fc,l)))e- E ^ Cl ^e^ )fc0 ^^ 

which gives us the equicontinuity. 

Arzela-Ascoli now gives us a subsequence of {/;} which converges uni- 
formly on compact sets to a function / : C — > C. We now show that this 
sequence can be picked such that / is actually a homeomorphism. To this 
end consider the inverse functions g\ = fjT . These functions satisfy the 
estimate: 

\gi(z) -gi(w)\ < 16tt 2 — 



log(e+^) 



Since 



i+lw(C) 



HMO 



25) almost surely we 



< K(C) and K E L 1 (D) (see Lemma 
immediately have that the sequence {gi} is equicontinuous. In addition, the 



integrability of distortion leads to the conclusion that / G W,' (for a proof 



of this last fact see [AIM09], theorem 20.9.4). 
The moduli 

and e~ CiN " 



The modulus of continuity is given by the relation between p n Pn +bn ^ Nn 

□ 



4 Vaguelets 
4.1 Construction 

Consider a wavelet basis {&j,i} of L 2 (M) with mother wavelet <5 : M. — > M. 
Following Donoho ( |D95j ) set 

<f,(t) := — [ e^&iuyul-Wdu 
27r Jr 

The vaguelet <j> is the half integral of $ and satisfies the following properties 
(we can choose q by choosing a suitable decay for <£): 

1^)1 <Ci(l + |t|)-(* +1 \ i £ K, 
4>(t)dt = 
\cp(t) - 0(s)\ < C 2 \t - s\ 
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Following Y. Meyer ( jM90| ) consider the periodized functions: 

oo 

*j(0) := a" 3 j>(2*(0-fc)) 

— oo 

The functions {1} U {9j(0 — k2~ 1 )} J) o<fc<2J form a periodic orthonormal 
wavelet basis of L 2 (T). 

We now introduce the periodic vaguelet. For a function on the torus with 
Fourier series f(9) ~ £ f(n)e 27Tind we have (-A)/ ~ £(27rn) 2 /(n)e 27rin0 so 
we may define the operator (— A)^ 1 / 4 by 

(-Ap/*/(n) := — l=/(n) 

Define now the periodic vaguelet ^(0) by 

${n) := J_ £(n) 

where is the periodic wavelet. 
Define 

fo(n) := -=^^-(^72^)(n) (12) 
V 27r|n| 

The vaguelets j are periodized versions of the <f>j i and have essentially the 
same properties. While the wavelets are a basis for L 2 (T), the vaguelets are 
a basis for the space H^ 2 (T) of functions / which have mean zero, and half 
of a derivative in L 2 (T). In other words, / ~ ^ n? iQ a n e 2n,lfl and the norm 
on the space is: 



II 2irin6\\ 

\\^a n e \\ H i/2- 
4.2 Properties 

In this section we present some properties of the vaguelets that will prove 
useful later. We start with: 




C \i> jyl {e)\de < C2-j 

Jo 
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We give a few estimates which we will need later: 



Lemma 6. For any q > and any 6 € J with | j| = 2 and J 
[2-^,2-^ + 1)]^ J 



\M0)\ < 



C a 



(1 + |2J'0-Z|)9 



(13) 
(14) 



The expression \ 2 3 9 — I \ should be understood modulo 2 J and can be replaced 

Remark 1. The quantity dls ^ J ^ i s the number of intervals of size 2~ 3 that 
separate J and J on the torus. 

Proof. The proof is a simple computation using the decay of the vaguelet 4>- 

J'H-^-/ + ^-)|<2. 7TT | W : r 

k 



(l + \2ie-l + 2ik\)i+ 1 ' 



(15) 



We may describe J as the interval [2~H, 2~ 3 (l + 1)]. Then 2 3 6 - I G [Z - 1, 1 + 
1 — /]. |Z — Z| equals the number of dyadic intervals (modulo 2 3 ) of length 2~ J 

separating J and J on the unit torus (quantity denoted by d%s ^i J - J " > ). This 
gives also the largest term in the series above. All the other terms decrease 
very fast and their sum is dominated by the first. 

An identical argument works for the derivative tpj(6). □ 

We also have the following: 

Lemma 7. For the family of vaguelets defined above the following relations 
hold: there is a constant Co such that for all dyadic I, all m 



|J|>2" 



< C V6> 



|J|<2-">,J(23J 



< Co yo £ i 



(16) 



(17) 



where 31 is the interval formed by the dyadic interval I and its left and right 
neighbors of the same size. 
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Proof. We recall that tfij are periodic vaguelets (defined on [0, 1]). We will 
first prove that 



E rt^ d9 

\J\>2~ m 



(m + l)ln2 



7T 



<C 



(18) 
(19) 



and then deal with the pointwise estimate. We will prove this inequality by 
reducing the computation to the wavelets 3>j on the line. 

In the following we will replace the notation ipj by ^>,y where J = 
[2~H, 2~i{l + 1)]. The periodic vaguelet ipjj was defined by: 



i>j,i{n) :-- 



II 



=* i (.-Z2-i)(n) 



(20) 



where \P was the periodic wavelet and fyj its refinement to level 2 3 . We 
also have 

^j(^T2-j){n) = 2-- ? '/ 2 e- 2W2 " J $(2-%) (21) 

where <3? is the mother wavelet on K. 
We then have 

m 2 J — 1 „i m IP — 1 m 



j=0 1=0 



j=0 1=0 n^O 



j=0 n^O 



2ir \n\ 



We remark that for m large enough: 



E^i^- m )i 2 = E 



1 



2vr \n 2- 



|$(n2- m )| 2 2- ,n 



1 

27 



I$(QI S 

lei 



(23) 

By ~ we mean equal up to a small error (and all such errors add up to at 
most a constant). We will prove that 



1 

2^ 



\Ht)\\ * ^ 2 
~Tti — d * = 



(24) 



This fact follows from the construction of the mother wavelet <&. We recall 
the construction procedure (see |M90j, chapter 3, or [BNB00J, chapter 7). 

One considers a function rrao(£) (also called filter) with the following 
properties: 
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mo is continuous and 2ir periodic. 

K(£)l 2 + K(£ + 7r)| 2 = i. 

mo(0) = 1 and m (£) + on [-f , f ] 



For any such filter one considers the father wavelet (aka scaling function) / 
given by 

/(£):= n m °(|) ( 25 ) 

and the mother wavelet will be given by the relation 

3>(2£) := e-«mo(e + tt)/(£) (26) 
We may now proceed with our argument: 

l*(£)| 2 f |m (e/2 + 7r)| 2 |/(e/2)| 2 /■ |m (e + vr)| 2 |/(O| 2 



J ~ir dC= j & di= ! 



\z\ 



-d£(27) 



We make the following observation (using the second property of the 
filter mo): 

We have used the fact that /(2£) = ?tco(£)/(£) (which follows from relation 



25). We make a change of variable in the second term to get: 



2,+1 \mi d( = r K«^)i 2 i/(oiy + r 2 iMP„ (3 „) 

Adding these terms for j > — n we get: 

l/gjp^ = r KK+^mei'^ (31) 

2~ n S 72-™ S 

The product II^eN 7710 ^) converges uniformly on bounded sets of E (see 
e.g. [BNB00]) and thus / is continuous. Since /(0) = 1 we get: 

/ c <%->-rn2 as n-)-oo (32) 

J2~ n S 
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which gives us ( 24 ) . 

To prove the pointwise estimate we use the (inverse) Fourier transform: 

Since e - 2 ™' 2 ~ J = for all n ^ and 2 j for n = we get: 

|VhE^^*p->*(-'(-*)) (34) 



by the computation in the previous part of the proof. 

The second inequality in the lemma is a consequence of the decay of the 
vaguelets (Lemma [6]). □ 



5 The Gaussian Free Field and vaguelets 

Heuristically, the Gaussian Free Field is a Gaussian "random variable" on 
an infinite dimensional space. A precise and correct definition is more sub- 
tle. We first give a few facts about (usual) Gaussian random variables and 
then extend the concept to infinite dimensional spaces. We follow the pre- 
sentation in |S07| . where Sheffield gives a good introduction to the GFF. 

Let (•, •) be an inner product on M. d and let \i be the probability mea- 
sure e~^ v ' v ^ 2 Z~ 1 dm(v), where m is Lebesgue measure on M d and Z is the 
normalizing constant. 

Proposition 1 QS07J). Let v be a Lebesgue measurable random variable on 
M. d with inner product (•,•)• The following are equivalent: 

a v has the (Gaussian) law \i. 

b v has the same law as Ylj=i a j v j where v\, . . . ,Vd are a deterministic 
orthonormal basis o/R d and otj are i.i. d. Gaussian random variables 
with mean zero and variance one. 
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c The characteristic function (Fourier transform) of v is given by 




d For each fixed w G M. d , the inner product (v, w) is a zero mean Gaus- 
sian random variable with variance (w,w). 

The Gaussian Free Field is supposed to be a variable on the infinite 
dimensional space H l {D) = Wq' 2 (D), where D is a subdomain of M. d . If D 
has no boundary the space = Wq' 2 (D) stands for the Sobolev space 

of functions with mean zero and one derivative in I? . This space is a Hilbert 
space with inner product (/, g)y = J VfVg. 

Ideally one would consider an orthonormal basis bj of this space and 
declare the GFF to be the random variable given by YlJLi a j^j with ctj ~ 
N(0, 1) i.i.d. However, this sum doesn't converge in H 1 and one has to 
consider its convergence in a bigger (Banach) space. 

Alternatively, one can define the GFF as being the formal sum h = 
Y^jLi a jbj with ctj ~ N(Q, 1) i.i.d. and {bj} an ordered orthonormal basis. 
For any fixed / € H 1 with f = Yl fijbj one can define the inner product as 
a random variable (h, /)y = lim^oo Ylj=i a i^r 

The Gaussian Free Field then becomes a collection of mean zero Gaus- 
sian random variables {(h, /)v}/g/f 1 with variance (/, /)v and covariance 
strucure given by 

%/)v(M)v] = (/,3)v 

On a manifold with no boundary D in M. d we also have (pi,/92)v = 
— (pi, — A/92)l2. This allows us to define the GFF as being a collection of 
mean zero Gaussian random variables {(h, p)} p e(-A)H 1 (D) with covariance 
structure given by 

E[(h,pi)(h,p 2 )] = / p 1 (x)G(x,y)p 2 (y)dxdy 

J DxD 

where G(x, y) is Green's function on D (the inverse of the laplacian operator 
on D). 

The Gaussian Free Field we work with is the trace on T of the 2- 
dimensional GFF. In stead of being a random variable on the space H 1 , 

1/2 

the trace of the 2-dimensional GFF is a random variable on the space H Q 
of mean zero and half a derivative in L? . Formally, this can be defined as 

Hl /2 (T) = {f\J fd0 = 0, (^) V2 / e L 2 (T)}. (37) 
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We recall that {1, \/2cos(27rn6>), v / 2sin(27rn6»)} is an orthonormal ba- 
sis of L 2 (T). This implies that { i — cos(2-7rnfl), J — sin(27rnl9)} form an 

\J 7r|n| y 7r|n| 

1 /2 

orthonormal basis of H (T). This allows Astala, Jones, Kupiainen and 
Saksman j A.1KS091) to define the trace on T of the 2-dim GFF as the the 
random distribution: 

x _ A n cos(2Trnt) + B n sin(27rnt) 



Jn 

n=l 



where A n , B n ~ A^(0, 1) are independent. 

We can also consider a wavelet basis {4>i} for L 2 (T). The image under 
(—A) -1 / 4 is a vaguelet basis {i^i} for Hq (T). Then we have that up to a 
probability preserving transformation the GFF can be rewritten as 

X = ^A I ^i 
i 

where Aj ~ N(0, 1). The factor y/n appears in this expression because it is 



missing in definition ( 38 ) . From this point onwards we will include it in the 
notation ipj. 

One can see the equivalence of the representations in the following way. 

1/2 

If we have two bases {Fj}, {Gfcjfor H (e.g. coming from bases {fj}, {gk} 
of L 2 ) then 

j j k k j 

Since Ej(Fj,G k ) 2 _ 1/4 = 1 for all k then E^G*)-!/^ ~ iV(0, 1) so 

AjFj = E B ^ B k e N{0, l) 

j k 

up to a measure preserving transformation. 

In |AJKS09] Astala, Jones, Kupiainen and Saksman used a white noise 
representation for the GFF. Gaussian white noise is a centered Gaussian 
process, indexed by sets of finite hyperbolic area measure in the upper half- 
plane and with covariance structure given by the hypebolic area measure of 
the intersection of sets. The trace of the Gaussian Free Field on T was then 
expressed as 

H(x) = W(x + H) , x £ T where (39) 
H = Ux, y) G HI - 1/2 < x < 1/2, y > 2 tan(|vrx|)} (40) 

7T 
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The geometry of the set H allowed Astala, Jones, Kupiainen, Salesman to 
decouple the variables on different scales in their main probabilistic estimate. 
While the vaguelets are slightly more complicated, their tails decay fast 
enough to allow us a similar decoupling. 



6 Random measure 

6.1 Construction of the measure 

We write the GFF as 

X = Y J diM&) (41) 
i 

where aj ~ iV(0, 1) and {V'/} are the vaguelets defined in section |4j scaled 
by the factor ^/tt as we pointed out in section [5] 

Take t k = t c - k~~< and n k ~ (k + 1 )7 e ic(fc+l) 3 ^ +1 > 7 . Define S := 

and 

S fe+ i(0) :=S k (9)+ Y, (aM0)-^$(P)),Vk>0. 

/:|/|6[2" n fc+l,2-( n fc + 1 )] 

(42) 

Here ai are independent centered Gaussian random variables of variance 
t k+1 for \L\ G [2 _nfc + 1 ,2-( n *+ 1 )]. The sequence {t k } is increasing to t c = 2. 
Define 

dv k := e Sk{6) d6 (43) 

F k ■= [ du k = [ e SkW d8. (44) 
J [0,1] J [0,1] 

It is easy to see {F k } is an L 1 martingale and hence it has an almost sure 
limit Fq . We want to prove this martingale is in the space L log L to ensure 
the F k — > Fq in L . This and Kolmogorov's zero-one law imply Fq is almost 
surely nonzero. The subcritical case, t k = t < 2 was studied by Kahane 
(|K85J) who proved that the martingale is in L p for some p = p(t) > 1. 

In the next result we will repeatedly use the equivalence of the LP norm of 
a martingale to the LP norm of its square function. Let {M k } with Mq = 
be an L p martingale for 1 < p < 00. Define the martingale differences 

Afc := M k — Mfc_i and set S k := (X)i=i^f ) • The latter is called the 
martingale square function and captures the LP behavior of the martingale 
(see |B55] ): 

^E[S p k ] < E[\M k f] < C p E[Sl\ (45) 
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The constant C p has order of magnitude ^= and is independent of the 
martingale. We will apply the right inequality repeatedly in the case when 
p G (1, 2] in which case the function x — > x p is subadditive and the inequality 
becomes 

k 

E[\M k n<C p E?£A?] (46) 

i=l 

Heuristically, one can interpret this as saying that the martingale differences 
behave as if they were independent. 

Theorem 8. Take t k = t c - jfe-T and n k ~ (k + l)7 e c*(fc+i) 3 ^ fc+1 > 7 A. Then 
the martingale F k satisfies E[F k log(l + F k )] < C . 

Corollary 9. The maximal function F* = sup fc \F k \ is in L 1 and hence the 
martingale {F k } is an H 1 -bounded martingale (H 1 denotes here the Hardy 
space). 

Proof of theorem [#| We begin by obtaining LP estimates on the martingale 
differences. 

F k+1 -F k = f 1 e s ^ - e s ^dd = f e ^W (e Ez(a^rW-^!W) _ ^ 
Jo Jo 

The exponent of the second term can be written as 

I ' \J\=2- (n k+ 1 ) ICJ ^ " ' |J|=2- n fc +1 

(47) 

(with the obvious definition). In fact, we may denote Aj, by Aj where 
j G {1, . . . , 2 nfc+1 }. Set A) = 0. We may now write: 



F k+1 -F k = /' 1 ^W(^ +1 ^-l)dff 
Jo 

l=l Jo 

i=i Jo 



XftdB 

i 
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where Y\ := (e 1 — 1) and are mutually independent and independent of {Xi} 
(stand for the other two terms). The independence follows from the fact that 
the normal variables which appear in the definition of A\ are associated with 
the dyadic intervals which are subsets of the dyadic interval of size 2~( nk+1 ^ 
which corresponds to I. 

We also have E[Y t ] = 0. This implies {£jLo Jo *i*Wo"* +1 is a mar " 
tingale with respect to increasing L. This implies (via the martingale square 
function) for p £ (1,2]: 

E[\F k - F k+l \P] < c p V E[\ f 1 XftdBlP] (48) 
l J o 

We estimate the term X\Y[dQ in the same way since it can itself be 
thought of as a martingale in the following way. Denote by J{%) the collection 
of dyadic intervals I C J such that |/| > |J|2~\ Set also J(— 1) = and 

' i=i /ej(i)\j(i-i) ^ 

We may now write: 

i Tl/g _)_ ^ — T^k 1 / + \ / + \ 

/ W<? = V / jf l(e S/ 6J (o(^H^--^*?W)_ e E/ 6J («-i)(^W-^W)) c «) 

= V / Xje EjeJ ( i - 1 )( a ^ l(fl) ""^^ (0) )( e i: ^ J (o\ J (i-i)(^^( fl )-- ! * 1 ^W) 
i=0 ^° 

= Z, XiZi ti _iTi ti dO 

We use the index Z in stead of J because we will sum by I later and J is the 
Z'th dyadic interval of length 2 _nfc_1 . 

The random variables 7} j are mutually independent (with respect to i) 
and are also independent of {X/Z; .}q _1 . They also have mean equal to zero. 
By the same argument as before: 

£[| / x,y,d0H < Cp V e[\ / x,z M _iT M d0H 

Jo i=0 Jo 

In J(i) \ J(z — 1) there are 2* dyadic intervals of length 2~ nk ~ l which we 
now index by %. The associated random variables are centered independent 
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Gaussians. So we may write yet again: 
21 rl 



Jo x=i 

■si 



2 l f l 
. vn 

X= 



and 



z* 1 21 f 1 

E[\ IiVWl^^E^I X^x^x^H 
Jo x=1 Jo 

Now |V X | = | e «'x^xW-^ x W _ i| < e \^r x (0)-^ x (0)\ laix _ 

t ^2 1 ipi x (d)\\ipi x (d)\ by the mean value theorem. We know that \ipi x \ < c 
and with the obvious notation we have \V X \ < e cWx W x \ipj x \. It is easy to 
see that E[e pcWx W x ] < C some universal constant. 
So we have 

E[\ f 1 X^^U^V^ef] < E[\ f 1 X l Z l ^ 1 U x - 1 e cW *W x \il> Ix (6)\dO\ p ] 
Jo Jo 



The last inequality holds by Jensen's inequality, while the equality holds 
because the random variables Xf, Zf i _ 1 , U^_ l: e pcWx W x are independent. 
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One may easily see that : 

[ |V/ x P<Ci|J x | 

Jo 

E[e pcW *WP] < C, universal constant 



E[Xf] = £[e pSfc ]e Ej=oE/CJ ; (p p) ~ 





2 


'fcH 


-l^?(fl) 




2 


'fen 


_l"/> 2 (9) 


'fen 


2 




2 


*fe- 





Putting all this information together we get: 

E[Xf]E[Z?._ 1 ]E[UV x _ 1 }E[^ w xW%] 

V fa 2 V i fo 2 si 'Hi^W 

< (J e ^\i\>2- n k VP P> 2 e ^|/|e[2- n fe- I ,2- n fc" 1 ]'- p P> 2 

for all /, i, x in their respective ranges. 
We apply inequality ( 16 ) to get 

e ( P 2 -p)^±iE |I|e[2 -n fc -i j2 -n fc -l ] ^(0) < Ce ( p 2_ p)4 l n2 ^±i 

Similarly, 

g^|/|>2-™fc^ W 2 < Q e (P ^-P) E m =l 2 



Finally we may write 



Jo 



k ( " a ~" a - l)ia " l2 e C (p 2 -p)»ln2^±l 
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Taking all this into consideration we get: 

2 n k+ 1 n k+1 -n k 2* 

E [\F k+1 - Ftf] < c 3 p CC p E E E 2-^+^2(*> 2 -^=i — "T l "- m 

1=0 i=l x=l 

< clCCf2 nk 2- pnk 2 ( - p2 - p ^°=i 2 E 2 i 2 _pi 2( pa - p ) i a 

i=l 

n k+l~ n k 

<c 3 CC p 2 -Mp-i)Hp 2 -p)i:Li ( ^r 1} ^ v 2- < ( p - 1 > + ( pa - p > i ^ 

i=i 

n k+l~ n k 



< c 3 CC p 2- nk( - p - 1 ^ 1 -2) V 2 -i ( p-1 )( 1 



i=i 

We have used the fact that 

E(n m — n m -i)t m _ n k t k v n m (t m — i m +i) ^fcifc 

m=l m=l 

which follows from the fact that the sequence {t k } is increasing. 

We make two observations. First, notice that t c = 2. Secondly, the first 
martingale difference, F\ — Fq, dominates all others and satisfies 



C 

2 ( pi -i)(in^i) _ 1 



Em - F n < - iw — (49) 



The power p\ is chosen such that (jpi — 1)(1 — Pj ^ 1 ) > 0. The closer t\ is to 
t c , the closer p\ is to zero and the worse this bound is. 

The next step is to obtain bounds on E[F k log(e + F k )\. 

By the mean value theorem (applied to f(x) = xln(e + x)) we may write 

E[\F k+1 ln(e + F k+1 ) - F k ln(e + F k )\] 
< E[\F k+1 - F k \ sup (ln(e + £) + 1)] 

< E[\F k+1 - F k \ ln(e + sup Fj] + E[\F k+1 - F k \] 

i<k+l 

We need to bound the first term. We start by applying Holder inequality, 
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and continue on the third step by Doob's inequality: 

E[\F k+l -F k \\n(e + sup FJ] < E[\F k+1 - F k \ p ] l l p E[He + sup F^} 1 /" 

i<k+l i<k+l 
poo 

E[ln(e + sup Fi) q ] = l + q A 9_1 P(ln(e + sup F<) > X)d\ 

i<k+l Jl i<k+l 

P(ln(e + sup Fi) > A) = P( sup F { > e x - e) < M^±ll ~ e -A 

i<fc+l i<k+l e — e 

/»oo 

£[ln(e + sup Fi) q ] = l+q X q ~ 1 e~ x d\ < Cq(q\) 

i<k+l Jl 

E[\F k+l -F k \\n(e+ sup F,)} < CE[\F k+l - Ftf] 1 /^.)) 1 '* 

i<k+l 

< CE[\F k+l - F k \P]^q = CE[\F k+l - F k \^-^. 

Now we put everything together (C is universal, i.e. it doesn't depend on 
k,p). 

E[F k+l ln(e + F k+l )\ < E[F k ln(e + F k )\ + E[\F k - F k+l \] + 

+CE[\F k+1 -F k \P} 1 /P-P- J 

We want to make sure that the second and third term add up to less than 
s = (fcjxp ■ We use the estimates above for p = p k +\, which we specify 
below. It depends on our choice t k . 

E[\F k - F fc +i| Pfe+1 ] < c Pfc+1 C2 nfc 2 (p fe+i" Pfc + l)E '=i 2 2~ nkPk + 1 = 

where we have denoted S k = Yli=i (jli ~ n J L ~ 1 " >ti . This implies: 

(E[\F k - F k+1 \ Pk + 1 ]) 1/pk+1 < C2 nk ^ k ~^ +{Pk+1 ~ 1)Sk+ ^ k ~^~ 
For t k = t c — k~ J we take p k = 1 + kr 1 '/2. It suffices to have n k such that: 

2 nfc ^Sl i+(pfc+1 " 1)5fc+!!: Sr 1 " 21n(pfc+1 " 1) <r — 

(A; + l) 2 

which is equivalent to: 

n k l ~ Pk+l + ( Pk+1 - l)S k + - 21n( Pfc+1 - 1) < -21n(* + 1) 

Pk+i Pk+i 
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A closer look at S k reveals that this is less than n k t k /2 so it suffices to find 
n k such that: 



n k h [pk+i ~ 1)— H — -21n(pfc+i - 1) < -21n(fc + 1) 

Pfc+i 2 p fc+ i 

t fc 1 \ , lnc Pfe+1 



n fc (p fc+1 - 1) f < ^ + 2 ln(p fc+1 - 1) - 2 ln(fc + 1) 

Which is equivalent to 



n k > 



ln ^+21n( K+1 -l)-21n(fc + l) 



(p fc+1 -l)(^--L 



2 Pfc+i ' 

lncp - 2pfc+i hi(pfc+i - 1) + 2p fc+ i ln(fc + 1) 

n k > 2 w s ■ 

(Pk+i ~ 1)(2 -Pk+itk) 

Keeping in mind that c Pk+1 ~ ^ 1)2 (the constant that gives comparabil- 
ity between the martingale and square function, squared) we see that our 
initial choice of n k satisfies this requirement. This concludes the proof of 
the L log L integrability of the martingale F k . □ 

6.2 Properties 

We now list the main properties of the random measure v = lim^ u k (see 



( 43 ) for the definition) . 

Theorem 10. The limiting measure satisfies: 

(a) Almost surely, for all intervals I, v{I) > 0. 

(b) Almost surely, for all intervals I we have v{I) < \I\ ak , where a k 
corresponds to the t k for which \I\ G [2~ nfc , 2~ nfc - 1 ). In particular 

„(/) < e ~\l^\. 

(c) For any subinterval I the random variable v(I) has all negative mo- 
ments. 

In particular, this measure is non-atomic and non-zero on any interval. 

Proof. We begin by arguing that the measure is a.s. non-zero on each inter- 
val. Consider an interval /. The event v (I) = is independent of the behav- 
ior of any finite number of levels in the GFF and hence a tail event. Since 
the martingale v k {T) is uniformly integrable we must have P{v k {I) = 0) = 
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by Kolmogorov's zero-one law. Thus the measure v is almost surely non-zero 
on the interval /. Since there is a countable number of dyadic intervals we 
get that almost surely the measure is non-zero on any such interval /. As 
any other interval contains a dyadic interval we get (a). 

We prove (b) by first showing that given a dyadic interval / with |/| = 
2~ n and n G [2~ nk , 2~ nk - 1 ) , then 

P(u(I) > \I\ ak ) < \I\ 1+€k . 

We have 

oo 

p{v{i) > |/|°*) < p(u k (i) > |/rv2) + J2 p (^+i0o - > i/r^*-*) 

i=k+l 

The terms in this sum can be bounded using the same estimates as in the 
proof of the uniform integrability. We obtain the following: 

P{Vk{I) > \I\ ak /2) < 2 Pk \I\ Pk ~ akPk 2( P k~ Pk ^ V:J 2 K = 4\I\<Pk- a kPk 

We use the notation ( p := p — (p 2 — p)|. The power pk corresponds to the 
maximum variance, t k , that appears in the definition of the measure v k . 
We can pick a& small enough to have Q Pk — a^Pk = 1 + £fc • 
The other terms satisfy the following inequalities (by the same estimates 
as above): 

P(\u i+1 (I) - Vi{I)\ > |J| a fc2 fe_i ) < C w 2 (i-fc)pi |J| Pi ~ 0fcPi 2 (ni ~ n)(1 " Pi) 2^- pi ^. 
We now check that for each i > k + 1 we have: 

C p .2^ i ~ k ^ p ^I\ Pi ~ akPi 2^ ni ~ n ^ 1 ~ Pi ^2^i~ Pi ^ 71 ^ L < 2 i ~ fc 2 _rifc ( 1+efc * > 
It suffices to have these inequalities for n = n^: 

C p .2^ i ~ k ^ Pi 2~ nk ^ Pi ~ akP ^2^ ni ~ nk ^ 1 ^ p ^2^~ p ^^ < 2 i ~ k 2~ nk ( 1+fLk ^ > 
This is equivalent to the following relation on rij {i > k + 1): 

(i - k)pi - n k (pi - a k Pi) + (n, - n k ){l - pi) + (pf - Pi)^ 1 <i-k- n k {\ + e k ) 
i(pi - 1) + n,(l - pj) + (pf - Pi)^7T - k< ^i - 1) - ^fe(efe + afePi) 
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For our choice of rii these relations are satisfied. To conclude the proof of 
part (b) we consider the following (with |/| = 2~ n ) 



rifc+i-l 



< OO 



1 - 2" £ fe 

If we shift the dyadic grid by 1 /3 we get good estimates for all intervals and 
then Borel-Cantelli implies (6). 

We now turn our attention to the existence of negative moments (part 
(c)). We will address this by estimates on the Laplace transform of the 
measure v in the vein of |AJKS09] . We start by analyzing f([0, 1]) and 
obtain estimates on the Laplace transform by finding a recurrence relation. 

Ideally, we would write 

K[0,1]) > B(M 1 + M 2 ) (50) 

where Mi and M 2 are independent and have the same law as z^([0, 1]), up 
to scaling. ^([1/4,3/8]) and ^([5/8,3/4]) would be candidates. 

Unfortunately, the two candidates are not independent and the law of 
^([1/4,3/8]) differs from that of z/([0,l]). The differences are: 

• The definition of ^([1/4,3/8]) contains vaguelets on levels {0,1,2} 
which are above (in dyadic tree sense) [1/4,3/8]. 

• The variances appearing in i/([l/4, 3/8]) are slightly different from the 
ones in the definition of f ([0, 1]) relative to the corresponding interval. 

• Modulo the first two differences, the law of ^([1/4,3/8]) differs from 
the law of ^([0, 1]), by a factor of 1/8 due to the change of variable 
that maps [1/4,3/8] to [0,1]. 

We start from the two candidates and decouple them in such a way that 
we get to our goal. 
First consider 

ee[0 'V&- 2 2 

Then ^([1/4,3/8]) > e x D([l/4, 3/8]) and i/([5/8, 3/4]) > e x P([5/8, 3/4]), 
where the tilde stands for the measures constructed starting only with dyadic 
intervals of size less than 2~ 3 . It is easy to see that X satisfies the same 
distributional inequality as the one given in lemma 12 Secondly, denote 



30 



I\ = [1/4,3/8] and I2 = [5/8,3/4] and consider the measure fj(ij) formed 
by using only the vaguelets starting with dyadic level 3 and corresponding 
to intervals J C 3I{. Define also 

Xil= ^. E MM-^AO))- 

e 1 J<£u it \J\<\h\ 

Then 

v{h) > e Xi Vi{Ii) 
The random variable Xj can be written as 

Xi = inf (Fi(e) - Fi(e )) + F^Oo) - sup V ^j(9). 

tieh J£3Ii,\J\<\h\ 

where F^Oq) = Yl\j\<£Mi a J' l Pj(^o)- This term is a Gaussian random variable 
with mean zero and variance X]|j|(23/- V'j(^o) which is bounded by a universal 
constant Co according to lemma [7j The last term is also bounded by Co , 
while the first is controlled by lemma [i~2j Then 

P{Xi > A) < c(l + \)e~ x2 / 2 

The random variables are now independent, but they are formed 

only using the vaguelets corresponding to 3/j and thus are quite different 
from ^([0, 1]). We make them more similar by the following argument. Let 

Y % := - sup V (bj^j(6) - %^ 2 j{9)). 

^ ■ J£3Ii,\J\<\Ii\ 

We take the random variables bjj to be N(0,tj) and independent of a j and 
of one another. Thus Y\ is independent of Y%. This allows us to write: 

Viiji) > e Yl v ifi {Ii) 

where z^o is defined by using all intervals J C [0, 1] and, when J <f_ 3/,, 
by using the random variables bj^. The two random variables ^1,0(^1) and 
^2,0(^2) arid independent. It is immediate that Yi behave exactly like Xf 

P(Xi >^)< c(l + X)e- X2/2 

Finally, we consider 

' k=l 1=0 |J|=2-"fc- i 
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where a^j ~ N(0,tk — ifc-i). The distributional properties of the field Zi 
are the same as for Xj and Yi by using the same argument in lemma 12 and 
observing that the variances form a telescoping series. 
We then have 

where 1^,1 are independent and have the same law as |^([0, 1]). 

All this work takes care of the differences we mentioned above. We are 
now in a position to write: 

M > e x e Xl e Yl e Zl -M 1 + e x e X2 e Y2 e Z2 -M 2 (51) 
8 8 

where M, Mi,M 2 have the law of ^([0, 1]) and the last two are independent 
of one another. 

Let C(s) = E[e~ sM ] be the Laplace transform of M. The recurrence 
relation ( |51[ ) implies 

2 ^ C^X^YuZi) + £2 
s 

where C(X, Xj, Yi, Zi) is a constant which depends on the corresponding 
variables. However, since all the variables X,Xi,Yi, Zi behave like Gaussians 
with mean zero and variance < 8, we won't worry about them. Setting 
f( s ) = ~T/2 + £( s ) we § e t f{s 2 ) < f 2 {s). We claim that there is an so such 
that /(s ) < 1 - e < 1. 

It suffices to find so such that £(so) 1S ^ ess than 1. 

We begin with some estimates (using notation M = -Foo) : 

P(M < 1/4) < P(F 1 < 1/2) + P(F 1 > 1/2, < 1/4) (53) 



The first term is bounded using (49) and lemma 11 



1 

p-i 



1 / 2 (p- 1 )( 1 -f)-1 > 
P(F 1 < 1/2) < 1 - — < 1 - — ] (51) 



where p is the value which corresponds to the variance t < t c appearing in 
the definition of F\. We denote the bound by 1 — 4e. 
The second term can be bounded by 

00 

P(F t > 1/2, Foo < 1/4) < J2P(\Fi-Fi^\ > 1/4*- 1 ) (55) 

i=2 
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and this in turn is small in comparison to P(F\ > 1/2) by the proof of 
theorem [8] and the choice of the {n^}. 

All this implies P(M < 1/4) < 1 — 3e. We now pick sq sufficiently large 

to have e _s<)//4 < e and < e. Then we have /(so) < 1 — e. We iterate 

V 

and get /(sq ) < (1 — e) 2 and by monotonicity f(s) < Cs~ s , where C < Sq 
and 5 ~ — — — — 



ln(s ) 
It follows that 



(56) 



We finish the proof as in [AJKS09 j by bootstrapping using the inequality 
between the geometric and arithmetic mean. We get 

E[M- 2 i] < C neg E[M-i} 2 (57) 

The constant C neg comes from the negative q moments of the variables 
X,Xi,Yi,Zi. It is a universal constant. 

To bound the first negative moment we need to do this operation ~ log 2 \ 
times. We get the bound: 

£[M- X ]< (C neg s 5 y = C^so (58) 

We give a summary of the variables that are relevant for this inequality: 

f (p-l)(l- P t/2)ln2 ^ 1 ln(l- C ) 

V C2P J ' e 2 ' lne 1 ' 

If we construct the measure v starting with variance t& = 2 — the bound 
reads: 

e~k-^ k \s ~ k^ k \5 — — (60) 

As we take tk — > t c this bound blows up. This points to the fact that as we 
approach criticality the measures become more and more concentrated. 

To obtain the negative moments of for some smaller interval one 
decouples the levels "above" /. These levels will form a centered Gaussian 
field which behaves like a Gaussian variable with variance ~ log 2 jjj and 
mean given by a sum of squares of vaguelets. For the decoupled part one 
applies the analysis here keeping in mind that the first variance has changed 
and hence the bounds have increased. □ 



33 



The next lemma appears in |L05j, but we give it here for the sake of 
completeness. 

Lemma 11. Let X be a positive random variable with E[X] = 1 and 
E[X P ] < oo for some p > 1. Then 

P(X > > (61) 

Proof. The proof is the same as in |L05j . It is immediate that E[X;x > 
1/2] > 1/2. Then 

E[X P ] > E[X p ;X>l/2] 

= P(X > 1/2)E[X P \X > 1/2] 
> P(X > l/2)E[X\X > l/2f 
E[X;X > l/2] p 



> P(X>l/2 

> 



P(X > 1/2)p 

1 



2PP{X > l/2)P~ 1 



□ 



Lemma 12. Let L be a dyadic interval of size 2~ l and consider the Gaus- 
sian field F = Y2j<^3i\j\<\i\ a J' l ( J j(^) with aj ~ N(0,tj). Then there exist 
universal constants c and C such that for fixed 9q 

P(sup \F(0) - F(Oo)\ > Cu) < c(l + u)e~ u2/2 . (62) 
eel 

Proof. The lemma is a consequance of the Borel-TIS inequality (see e.g. 
|AJKS09| ) for the field F(6)-F(6 ). We merely need to prove that E[\F(6)- 
F{e')\ 2 ] < L\6 - 9'\ for 9,6' C /. To this end observe that (0* G [9,6'} and 
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depends on J): 

oo 

E[\F{9) - F(6')\ 2 } < E tj \ 9 ~ 0'\Wj{e*)\{\M0)\ + l^(^)l)(63) 

j=i \J\=2-i ,J<£3I 

<c\e-e\J2 E (i + we* - mi- 1 (i + we - W q (64) 

3 \J\=2~i ,J£3I V 1 17 

<C|fl-0'|V2 J V , . , * - rXIN9 : (65) 

<c\o-*\Y,* E (i+^^^-^lEa?^^^) 

< C\e-ff\2 i {(a't) 

since we have g > 1. So L = C2\ Since we are interested in the supremum 
over an interval of size |/| = 2~*, Borel-TIS inequality gives the desired 
estimate. □ 



7 Decoupling 

In this chapter we prove that there are sequences p n , p n , N n , b n , c n such that: 

PfcMod(G(A(z,p n pl2p n pti)) < c n N n ) < p nP £ +K)N " 
i=i 

for any z£T and any mapping G with Beltrami coefficient p. 

Before we proceed we recall a few notations: in the following we will work 
with dyadic intervals /, J. V n is the collection of dyadic intervals of size 2~ n . 
For a dyadic interval /, is the union of / and it's two neighbors of the 
same size. We will extensively use the notation Cj := {(x, y)\x G /, 2~ n_1 < 
y < 2~ n }. These sets are called Whitney squares/boxes. 

We now fix a point z G T and n. For any i G 1, . . . , N n define annuli 
Ai := A(z,p n p l n ,2p n p l n ), and balls B { := B(z,p n p l n ), B\ := B(z, p n p l n / '4) . 
For simplicity of exposition, we will take Ai to be square annuli. Because 
the picture is a local one, we think of the annuli as being centered at a point 
on R. 

The annulus Ai can be divided in two parts: 

• Ri t i = Aid (U/gCjC/) j where Ci is the collection of all dyadic intervals 
/ of size at most p n p x n j\ such that / D Ai ^ 0. 
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• Ri,2 — Ai \ Ri : i 

A big modulus mod(G(Ai)) is a consequence of controlled distortion in 

Distortion in any Cj depends on the doubling properties of the random 
measure v in Let / be a dyadic interval whose doubling properties 

affect the distortion in A{ and any J = { Ji, J2} with Jj C 

For k G { — 1, 0, . . . ,i — 1} define 



k,k :=sup ^ ajipj(9)-M ^ ajipj(9) 

8Bi J.CjCWB, 
By lemma [7] we have 



8Bl J:C7jCl6B fe ,Cj^l6B fe+1 " J:C*jCl6 J B fc ,Cjfz:i6B fc+1 



A- 



l/(J 2 ) fj(J2) 



where Vi is the measure obtained only by using the dyadic intervals which 
are included in 16-Bj and C is a universal constant coming from the terms 

The t, satisfy the following: 



Lemma 13. For any indices i, i' the random variables tj & and ty &/ are 
independent if k ^ k' , and satisfy 



P(ti,k > Cuy/pUpn (k+1) = Cupt h - 1)/2 ) < c(l+n)e-" 2 / c , k G {-1, 0, . . . , 

where c, C are universal constants 

Proof. Consider the Gaussian random field: 

X{6) := Yl a ^/(#) " Yl a lM8o) 

\I\=lGp nP k n +1 \I\ = 16p nP k n +1 

By the mean value theorem E[(X(9)-X(8')) 2 } < Cij|0-0'| £^_ nP 'l fc+i l^r 1 (lV'/(^)l + 

\I\ — lQpnPn 

\ipi(8')\). Using the properties of the vaguelets this is less than Cti\6 — 
®'\Pn l P~k+v We look at this random variable on the interval R n 8Bi which 
has length 8p n p l n . By applying the Borell-TIS inequality we get the desired 
distributional inequality for supX(#). 

The same inequality holds for the sup of the field where we use —aj in 
stead of a/. The 9q terms cancel and we get the conclusion. □ 
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7.1 Distortion in Rii 

The region Rii has two pieces R il and R\i- We describe the decoupling 
for the left piece; the decoupling for the right piece is analogous. 

Let C\ i be the colection of all dyadic intervals I which do not appear in 
any of the to-, I C WBiDM and which satisfy one of the following conditions: 

• |J| > p n pi/4 

• if \I\ < p n p l n /A then IC\Bi + \ = and / does not intersect the p n p^/4- 
neighborbood of R\ 1 . 

Let trli := sup(J2 j€C l a jV>j(#)) ~ ' mi (J2j e c l a J^j(°))^ where the sup 
and inf are over 9 € R\ i D R. These variables have the following property: 

Lemma 14. The random variables tr^ and try y are independent if i ^ i' , 
and satisfy 

P(tr i:i > Cu) < c(l + u)e~" 2/c , 
where c, C are universal constants 

For each k > i let Cj & be the collection of dyadic intervals of length 
at most p n Pn (on M.) in \ B^+i to which we add the dyadic intervals 
J that are inside B^ and intersect -Bfc+i, but are not contained in it. Let 
tn t k ■= sup(J] J6C . k ajipj{0)) - inf (Y^Jec ltk ajipj(0)), where the sup and inf 
are over £ R%,\- We may write then 

for each J = {J\, J2} which influences the distortion in the region R\ l . The 
random measure UiQ is obtained by considering only the random variables 
a j corresponding to the J which are subsets of the p n p l n /4- neighborhood 

of r\ a n m. 

The terms of the form ^2jtjipj(9) cancel in the quotient ^jj^y because 
of lemma 

We should actually consider the two parts of Ri t \ separately in the com- 
putations that follow. However, their behavior is very similar so we allow 
ourselves to treat them ftS ft unit. See figure [IJ 

The random variables trik are supremums of Gaussian fields and they 
behave like Gaussian variables. 
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Lemma 15. The random variables ir^fc andtfi^y are independent ifk ^ k' , 
and satisfy 

P{tr itk > Cup^~ i){q - l) ) < c(l + u)e- u ' 2 / c , k > i 
where c, C are universal constants 




Figure 1: The decoupled variables used to control distortion in regions 
Ri ; i (shaded) 



Proof. The proof is very similar to the one of lemma 12 
Consider the Gaussian random field: 



X(9) := "iMO) ~ E a ^/(#o) 



Set / = E- j fl R. For any interval J 6 C^k we have dist(I, J) > p n p l n /2. 
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We may now write {6* G [8, 9'\ and depends on J): 

oo 

e[\f(6) - F(e')\ 2 } < E tj\e-e'u'j(en\(\M0)\ + \M0')\)m 

i=-iog 2 pnp\ \J\=2-3,JeC i:k 

<C\9-0\Y^ (1 + \2i0* - ZD"" 1 (1 + \2i0 - ^ 

j \J\=2 'jJeCj,* 

j \j\=2-i,JeC iik y 1 v ' >u 

< c\e - e'lYv Y t — 1 k^ < c\o-9'\ V / . 2J N9 . (72) 

3 d=23p n p* n /4 J=-log 2 PnP„ 



(o n k ) 2q - 2 

< c\e-e'\ ) pnPn { . (73) 



since we have g > 1. Since we are interested in the supremum over an 
interval of size ~ p n p l n ^ Borel-TIS inequality gives the desired, estimate. 

□ 

The random measures z/^o are independent of one another. While we can 
not say that with high probability the distortion in the regions R\ 1 and R\ l 
is bounded (this was the case in [AJKS09] ). we will construct a stopping 
time region inside each of these where the distortion grows in a controlled 
fashion. 

7.2 Distortion in R it2 

Let Zj be the set of I 6 T> such that Cj intersects Ai and |7| > p n p z n - The 
distortion in Ri% is the same as the distortion in all Cj for / £ Zj. This is a 
finite (and universally bounded) number of intervals. For each of them we 
only need to control ~ 2 10 pairs J. We have already decoupled the influence 
of variables aj with J ouside of 16-Bj on the measure v. For each pair J we 
can write (following |AJKS09| and recalling that B\ = B(z, p n p l n /4)) 

^ (J) - ^ (J 2 \ B if J 2 \ B,) + <MBd + (74) 

ui(jj n sj) = 2 ^(Jj n \ 4) (75) 

fe=i+l 



39 



Define 



Li,i= Yl ^(Ji\^',J 2 \^)(76) 
^ UitJt n(BL , \ B')) 

U, k = V 1 rrv^k + (i <+2) fori + i < fc < ^(77) 

An upper bound on X^fe>j ^«,A; m eans the distortion in region Rip is 
bounded. Before we give a distributional inequality for these random vari- 
ables, we need to decouple them one more time. 

n{Ji^{B' k _ x \B' k )) ^, fc (Ji n {B' k _ x \ B' k )) ^_ mT rm) 

where measure z/j & is constructed using only the random variables aj for 
which J C l6Bi and J <f_ \B' k - In addition, the sup and inf are considered 
over the set E n (QBi \ B' k ). Denote these variables by Ylj>k s^kj, where 
S{j t k involves only the dyadic intervals J which are subsets of ^Bj_ 1 , but 
are not subsets of \B'-. 

For simplicity of notation, we will use for the sums above, but in- 
volving fjfc. Distortion will then be bounded by 

; > ( 79 ) 

k>i 

For a picture of the decoupled variables see Figure [2j 

The random variables Sj &j have similar properties as tr^j. above due 
to the same reason: they are supremums of Gaussian fields. 

Lemma 16. The random variables Sijy and Si',fe'j' are independent if j / 
f, and satisfy 

P(s tAj > up^ 12 ) < Ce- u2 '\ j>k>i 
where C is a universal constant. 



Proof. The usual argument involving the Borel-TIS lemma (see e.g. |AJKS 09j) 
doesn't give a good estimate, because the lemma deals with more general 
Gaussian fields. We will give the simplest possible argument. 

For simplicity of notation, let = p n p\ and 2~ l ° = p n pn. All the sums 
refer to intervals J which appear in a* & and 9 is any point in Mfl (6Bi\B' k ). 
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2B; 



16 B: 



Figure 2: The decoupled variables L^. control the distortion in the shaded 
region using the vaguelets from outside the green square. Sj^j control dis- 
tortion in Ri t 2 using vaguelets inside the two red squares. 



Then 

oo , 

P(J>jVj(0) > A) < £p( E I^H^WI ^ ( 8 °) 

J l=l \J\=2~> 

oo . 

<E E p (\aj\\M0)\ > ^^z^) (81) 

«=«0 \J\=2~ l 

°° (7 A 

^E E ^i(^^^^y) ( 82 ) 

°^ _ A 2 2-?('~'l) ~ A 2 2 ( g -2)(i-i ) + g(i t ,- il ) 

<<?E E e 4 22(i "' o) <<?E 2 e 4 (83) 

Z=Z \J\=2~' l=lo 

Here C is a universal constant and dist is the distance between the region 
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where J are and the region where we take the sup/inf over. The statements 
above should be understood for all 6 € Rn(6Bi\B' k ), and not for a particular 

e. 

Replace A by A2-^°-' 1 )/ 2 and get : 

P(sup V aji>j{e) > A2-^°-^/ 2 ) < Ce~ x2/i (84) 
9 J 

Recalling that = p n Pn and 2~ l ° = p n p J n we get the conclusion. 

□ 

The random variables satisfy the following distributional inequality. 

Lemma 17. There exists a n > (depending only on the sequence of vari- 
ances in the GFF) and C n < oo (independent ofi,k,p) such that 

P(L hk > A) < C n A-Vi fe - i - 1)(1+a " ) 

In addition, and Ljj are independent if k < j or I < i. 

Proof. We recall that Lj & is a sum of a finite (and universaly bounded) 
number of terms of the form ~^^~p^g?y~^ • F° r simplicity of exposition 
we will redenote this quantity by where /, J are two intervals of sizes 

1^1 ~ PnPn and | J| ~ PnPn' 1 - 

The measure v is constructed using only the random variables aj for 
which J C 16-Bj and J ^ ^-BjL- The distributional properties of this measure 
are no different from the properties of the measure constructed using all the 
J C 16-Bj. This is the case because the vaguelets corresponding to \B' k form 
a Gaussian field with a controlled variance when evaluated inside 165, \ B' k . 
We use the same notation, v, for this more "complete" measure. 

We now have 

^>X] < P^>A/2Wr*< J >-?( J >>A/2 



Here the variables are the martingale approximations to v. We have the 
following trivial inequality: 
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Let P2 > q 2 for which E[u 2 (J) P2 ] < oo. Inspection of theorem [8] reveals that: 
E[v 2 (J) P2 ] < C P2 \jMl\ (Az ^ (86) 



where C P2 :=p 2 - ( ^p^. 

We know that v{I) has negative moments of all orders. In particular, 
theorem 10 gives us the following bound for the negative moment: 

/ J_ \ 3292 

E[ v {I)-i^\ < C\I\~™ 2 C*l gSl (87) 



where q 2 = P2 lq 2 (the conjugate of ||), and <5i,si are the constants which 



appear in the proof of theorem 10 



Applying Hoelder inequality and combining the last few estimates we 

get: 

fvo(J) \ n f2\^ 2 - /■ I 92 (P 2 ~P2)*2 . / 1 N<?2 

P > A/2J < CC 2 p2 (jj |J|«Ww|j|« 2 |j|- 92 (2h Sl j (88) 

<c(a, ft) ti)M (89) 

<C(h,P2,tx)^0^) +O2 (90) 
where C(q 2 ,p 2 ,ti) contains all the constants and q 2 Cp 2 /P2 = I + a 2 . This 

, 92 



constant is dominated by y2 s i sij and 8\, s\ depend on the first variance 
which appears in the definition of v (hence the parameter t\). 
We deal with the terms of the form 

P \ W) > A/2 

in the analoguous fashion. To bound the numerator one chooses Pi+i for 
which the corresponding moment of Vi+i( J) — Vi{ J) exists. The computation 
of this moment is basically given in theorem [8| 

E[\u i+1 (J) -Vi{J)\ n+1 ] < C i+1 \J\ c n+x \I\ {Pi+1 ~ P 2 i+l)H 2 -K-i)(Pi+i- 1 )(l-P»+i^/ 2 ) (9 

where C i+l = 2 ( P<+1 -i)(i- P<+lt<+1 / 2) and |J| = 2~K 

Combined with the negative moment estimate this gives: 

p r ^-w-^ > V2 , + a < Cfe+1 , B+1 , tl) W 



!/(/) ' / vii-t-i'^-t-ii ^ ^ y |j|9 i+ iC Pi+1 M+i 
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So we get (after we change the index) 



1=2 1 1 1=2 I I 



Since |/| ~ p n p' l n and |J| ~ p n p k 1 we have 14 ~ p k n 1 * > 1 4 > 
p N ", 

rn 

The variance ti corresponds to levels ~ p n Pn n which in turn is larger 
than 2 _ni in the notation of this lemma. 

|/| is between ~ p n Pn n an d Pn- To get the desired estimate it suffices to 
have: 

C(g 3 ,P3,ti)2^2-" 2( ^- 1 )| < l - p *n ( M2 Ar n) feC P3 /P3-i) (gg) 

Since 2~ ni+n ° < p^ n (hq corresponds to the levels above p n ) and 2 - ™ 1 < 
p n Pn n i it suffices to take: 

C(g 3 ,p3,ti)2 3 *2" n2Kp3 " 1) S < 2 (-"i+"0)a2-n 1 a 3 (g 4 ) 

If we also have the following relations on n k : 

C(q k+2 ,Pk+2,tk)2 3llk+2 2 fc+lUpfc + 2 W < 2(-^+^-iK+i-^+2,VA: > 1(95) 



we get the desired conclusion. This is because these relations, although not 
the same as the ones present in the sum, dominate the latter. 

In the case when t k = 2—k^ &ndp k = l+/c _7 /2 we get £ Pfe = l+/c~ 37 /8. 
We take q k = p k (l - fc~ 3 T/16) and get a k = k~^ /IQ - k'^/U ~ k-^/16. 

Since C(q k +2,Pk+2,t k ) < C!f e ~" g k k l2lkl relation (95) can be simplified as 
follows: 

ck^ kl - n k+1 a k+2 < {-n k + n fe _i)a fe+1 - n k a k+2 (96) 

so it suffices to take n k+ \a k+ 2 ~ c(k + l) 3 T( fc+1 ) 7 or n k+ i ~ (k + 2) 37 c(A; + 
1 \3 7 (fe+ir_ 

□ 

Recall that we are trying to prove an estimate on moduli of annuli at 
scales between p n Pn n an d Pn- All the work we have done decoupling the 
distortion was done to address a fixed n. The distributional inequality is 
satisfied by all Lj k , where i < k < N n . The constant a n depends on first 
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variances in the definition on the measures Ui^. For each n, all the first 
variances are equal to the same value t n . As n increases, a n — > 0. This is 
a big difference between the critical case and the non critical case treated 
in [AJKS09J. In the sub-critical case the variables Lj & were defined for all 
scales k > i and they all shared the same value of the constant a. 

8 Random tree 

We need to control the distortion in the regions Ri2- The authors of 
|AJKS 09j were able to get bounded distortion with high probability. In 
our case this is not possible anymore essentially because as the variances 
tn — > t c we lose control over Li±- We use a stopping time algorithm to 
construct a random tree. We devise a collection of rules which we apply to 
dyadic squares. If all these rules are satisfied for a particular interval/square, 
then the distortion will be controlled. Our goal is to obtain an infinite d-ary 
surviving tree where the distortion is controlled. 

8.1 The rules that define the tree 

In the following A, At, B are large constants, 5 a small constant and N is a 
large positive integer. 

We construct stopping rules on a tree in which each node represents an 
interval of size 2~ l ( N ~ A \ For each interval /, a node in this tree, we denote 
by Ans n {I) the ancestor n levels above (and of size 2~( 4 ~ n )( Ar_4 )). 

We start by considering a dyadic interval / of size 2~ J ( Ar ~ 4 ) and measure 
v which is constructed only using the vaguelets corresponding to intervals J 
which are subsets of / and its closest four neighbors (two to the left, two to 
the right - call this set J {I)). Out of all the dyadic subintervals of of 
size \I\2~ N i we select and mark half in an alternating fashion and call them 

1 k- 

The interval / survives if all of the following good events take place: 
1) 

^\ll\<v{I l k )<A t \ll\,\lk (97) 

where t represents the variance (in the definition of the GFF) for which 
all levels with that variance are below the level of I° k . 
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2) For each J e {/, Ii,I r } ( 1 and r stand for left and right neighbours) 
J2v(li)<A\J\2-3 s ,Vj>l (98) 



3c ' 



3) For each n < i: 

sup V ojVjW - aj^j(^) < 52- n (99) 

06.7(1) Jc7(An Sn ((J))\7(/) 

inf £ oj^j(fl) - ajVj(^) > -B2~ n (100) 

eeJ(/) Jc7(Ans„((/))\.7(/) 

where 9j is the center of /. 

4) 

sup ajij>j{0)-ajil}j{Oi)<B (101) 

9ej{1) JCj(I),\J\>\ll\2S 

e ™L a.j1>j{p)-a J 1>j{p I )>-B (102) 

where 9j is the center of /. 

All the constants are chosen such that 

P(rules hold) rj 1 (103) 

The constants At vary with the level at which we apply the rules. As we 
apply these rules to smaller and smaller intervals /, the properties of the 
measure v become weaker and weaker. We want to make sure that the 



probability (103) doesn't change as we go deeper and deeper. 

If all these rules are satisfied the distortion inside between heights 
\ll\2 5 and 2~ n = 2~ i ( Ar " 4 ) = |/| is bounded by D n . We reiterate the fact 
that this sequence won't be bounded. 

If all the rules are satisfied for J, we then look at its "children". We 
consider all the dyadic intervals of size 1 1\ \ 2 4 inside the region j (I) . We don't 
have control over the distortion in the boxes corresponding to these intervals. 
We select one third of these intervals in such a way that any two selected 
intervals are separated by four boxes which we do not select. See Figure 
[3J There will be a total of 2 Ar_4 /5 such intervals. We now run the rules 
for each one of these intervals. When doing this, we consider the measure 
v constructed only using the vaguelets corresponding to J (this interval). As 
we run these rules for smaller and smaller intervals we obtain a random tree. 
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Figure 3: The interval in the middle is /. The black intervals are the marked 
The green boxes represent the area where the distortion is under control. 
At the next step we only look at the white boxes that lie between four red 
boxes(here shown as two due to lack of space). 



8.2 A surviving d-ary subtree with high probability 

The goal of this section is to provide estimates on the probability that each 
rule fails and then to obtain an estimate on the probability that there is a 
d-ary surviving subtree. 
We have 

Lemma 18. 

2 

P(rule 1 fails) < — 

Proof. It is immediate that 



(104) 



For the other inequality we need to use the negative moment estimate ( 58 ) . 
We first write as e^--|/i|£'([0, 1]) where e^* contains the vaguelets from 

levels above that of l\ and below that of /. Strictly speaking, v is slightly 
different from the original construction, but there is no difference relevant 
for out computation. 
This implies 

f(j;|lil>>l'i))< y ^ (105) 
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We recall that the constants satisfy: 



i 



/ (p-l)(l-pt/2)ln2 ^ Ine 
V C2P y ' ° e 2'° ln(l-e) 1 j 

where t = tk for which all the re^ levels with t = tk are completely below 
level 2~ l ( N ~^ +N . p is the appropriate power. 

If we set A th = AClleg So,k eN we get the conclusion. 

□ 

Lemma 19. 

C 

P(rule 2 fails) < — 
Proof. We start by obtaining an estimate for a particular j. We remark that 



and that the latter is simply the limit of a martingale G m = v m (Ujj jlfy 
constructed in the same way as the original Ff.. There are only a few differ- 
ences: 

• The measure v is constructed only using vaguelets ip* corresponding 
to levels starting at 2~ l = \ J\ and inside 

• The index m = stands for oIvq = d9, while v\ stands for the measure 
obtained using all Gaussian random variables and vaguelets having 
the variance of level 2~W JN I ~ i ) + o N ) (the level corresponding to I J k ), 
denoted in this proof by t%. 

We will obtain estimates on the following 

P(|G m+1 - G m \ > 6A|J J 2 2 3 * ) (107) 

A look at the argument in theorem [8] reveals that the same argument works 
in this case. The only differences come from the fact that we have to replace 
the interval [0, 1] by S = U^ cJ I^ (and the fact that the vaguelets decay fast 
enough) . 

We get the estimate: 

E[\G m +i -G m \ p ] < cC2- ip 2( N -V j 2- Njp 2& 2 - p ^2- nm b-W-^)(iQ 8 ) 
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The n m , t m are all relative to the dyadic level of S. 

We emphasize that the p corresponds to the t m , so the estimate should 
read: 

E[\G m+1 - G m \ Pm ] < C2" iPm 2 ( ^ 1)j 2"^ Pm 2 (p -- pm) ^2~™ m(Pm " 1)(1_a ^ m ^109) 
This immediately leads us to: 

P{\G m+1 - G m \ > 6AlJ j 2 / ) ^ n °) 

171 TT 

= C2~ i(1 ~ ipm) 2~^ (Pm ~ 1)+(p ™~ Pm) ^2 _nm( ^^ 
This immediately implies: 

P( v{l{) > A\ J\2~i s ) < C2-i( 1 - 5p ^A- 1 (113) 

which in turn gives us 

P(rule 2 fails) < OA' 1 (114) 

□ 

Lemma 20. 

P(rule 4 fails) < C(l + CB)e~ CB2/2 

Proof. Rule 4 is a condition about a centered Gaussian Field with covariance 
bounded by C\P k \~ l . We are interested in its supremum/infimum over an 
interval of size 3|7|, so by Borel-TIS these behave like Gaussian variables 

with mean zero and standard deviation C\J\P k \~ 1 \I\ = C2 N / 2 = C. 
This gives us: 

P(rule 4 fails) < C(l + CB)e~ CB2/2 (115) 

□ 

We remark that rules 1,2 and 4 alone (not considering rule 3) give rise 
to an independent tree. The probability that a node fails is given by pf = 

C +C ( 1 + CB)e -CBV2 

Lemma 21. 

P(there is d-ary subtree, which survives rules 1,2,4) > 1 — 2pf 
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Proof. The process under scrutiny is dominated by a Galton Watson process 
because the death rules for one child is independent of the death rules for 
another child and the death probabilities are uniformly bounded by the 
estimates we have. The maximum number of descendats is M = 2^ N ^ 4 ^ 5 . 

Let t be the probability that the GW tree T has a d-ary subtree. Let 
T m be the tree T truncated at level m (i.e. after m steps). Let q m be the 
probability that T m doesn't contain a d-ary subtree. Then 

r = lim(l - q m ) 

In addition we have the recurrence relation: 

q m = Gd(q m -i) 

where Gd(s) is the probability that at most d-1 children are marked, when 
we mark each child independently (of each other and of the initial GW birth 
rules) with probability 1 — s. 

We have the following bound on Gd(s): 

d-l 

G d {s)<Pf + Y,C ] M s M ^ (116) 

3=0 

If the probability, pf, that an individual dies is small to start with (on the 
order of C^ /2 {2p f ) d ~ 2 < 1), then 1 - r < 2p f . □ 

Another essential estimate is the following: 

Lemma 22. P(there is a d-ary subtree, which survives rule 3) > 1 — ps 

Proof. We divide rule 3 in infinitely many rules. Rule n is: / survives if 

sup Yl aj^j(9) - aj^j) < B2~ n (117) 

eej(i) Jc j(Ans n ((i))\j(i) 

sup J2 aj^ J {e)-aj^j{e I )>-B2- n (118) 

eej(i) Jc j(Ans n {(i))\j(i) 

where 9j is the center of /. Notice that for the intervals of generation 
0, 1, . . . , n — 1 the rules are satisfied by default. In addition, the intervals on 
level n are pretty much perfectly correlated with one another. The intervals 
on level n + 1 are correlated in groups of 2^- 4 )/5, but these groups are 
independent of one another. 
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Let T m be the tree after m steps. Let q n ^ m be the probability that T m 
has no d-ary subtree which survives rule n. We seek a recurrence relation 
on q n>m . 

If T m has no d-ary subtree, one of two things must have gone wrong: 

• the descendants up to level n die (call this event Bad\), or 

• the descendants up to level n survive, but less than d of the generation 
1 descendants have a d-ary subtree (call this event -Baa^). 

There are a total of M = 2 N ~ 4 /5 possible descendants in generation 1 
and M n in generation n. The probability one of these dies is 

B 2 2 2n((AT-4)/2-l) 

P(sup(Y,) > B2- n ) < (1 + B2 ni - {N ~^/ 2 -^)e 3 = : Vn (119) 



This follows from an argument similar to the one in lemma (13). It's a 
consequence of the Borel-TIS inequality applied to a centered Gaussian field 
which covers dyadic levels up to ~ 2~( n_1 )( Ar_4 ) over a set of size ~ 2~( N ~^ 
(and behaves thus as a centered Gaussian random variable with variance 

C2 -n(N-4)/2y 

We can now conclude that P{Bad\) < p n M n . 
The second bad event can be bounded as follows: 

d-1 

P(Bad 2 ) < J2 Cnten.m-l)^. (120) 
i=0 

This holds for two reasons. First, once we condition on the survival of gener- 
ations 1, . . . , n, the subtrees starting at generation 1 nodes are independent. 
Secondly, although they don't have exactly the same law as the original 
tree, they are stochastically dominated by a tree which follows the same 
rules with random variables aj ~ N(0,t c )- 
We may thus write 



R2 9 2n((JV-4)/2-l) / oiV-4 \ n 

< (l + B2^ N -^)e- B2 " 2 (2_J +£c;(g n , m _ 1 )< 



\N— 4\ n d— 1 

<ln,', 

\ : I / 

i=0 

(121) 

One can easily see that this recurrence relation implies that q n = lim m q n ^ m < 
2p n M n , which in turn implies: 



P(no d-ary subtree, which survives rule 3) < 2^2 Pn M n (122) 

n 

which can be made arbitrarily small by choosing B large enough. □ 
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Putting the last couple of lemmata together and using the notation p = 
P3 +Pf we get: 

Lemma 23. P(there is a good d-ary surviving subtree) > 1 — p. 

Notice that if we take A ~ e 5B the probability p < §. By "good" we 
mean that the distortion is controlled. Explicitely, the distortion is Dt k+1 < 
CAA tk+1 2 N e 5B ~ CA 2 A tk+1 for levels [2" nfc , 2~ n ^ +l ]. C is a universal 
constant (depends only on the properties of the vaguelets). 

9 Modulus estimate 

In this section we are concerned with a deterministic modulus estimate. This 
modulus estimate is similar to some of the estimates of J. C. Yoccoz on the 
local connectivity of the Mandelbrot set (see e.g. [H92], [MOOJ). 

Assume we have a mapping G : C — > C with distortion \i and an annulus 
A centered on the real axis (for simplicity we will use a dyadic annulus of 
radii 1 and 2). Assume fj, = in the lower half plane. 

Assume the distortion in the top part of the annulus is bounded by D 
and that in the two sides we have connected sub-domains fij. and of A 
flj is constructed in a stopping time fashion on a dyadic grid. At level 2~ N% 
the number of surviving intervals is d l and the distortion there is bounded 
by D{. One should think of Qj converging to a Cantor set Ej on M Pi A. See 
figure [4j 

The next theorem is a generalization of the following statement: if G is 
conformal in the top of A and inside Q\ U O2 then mod{F{A)) > cq which 
depends on the logarithmic capacity of the sets E\,Ei- 

Theorem 24. Mod(G(A)) > a = const. 

Proof. One way to prove this result is to construct two closed curves inside 
G(A) which are neither too long, nor too close to one another. When G 
is conformal (as above) this is a consequence of Pfluger's theorem (see for 
example [BB09J for a reference) and a modulus estimate. 

We will use a more direct argument here. We know that modG(A) = 
modG(f) ' wnere r is the family of curves which connect the two components 

of the complement of G{A). We want to obtain an upper bound on modG(T) 
and we do this by constructing a good metric in G(A). 

We consider the metric p on G(A): p(w) = r^r^jn where G{z) = w 
and u is defined below. 
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GittOfiPOririal how 



Figure 4: The annulus A and the mapping G. The black regions are the 
ones where we have no control on the distortion. 

Then modG(T) < J J n p 2 dA(w) < J J G(n) \Vu\ 2 D(z)dA(z), where D(z) 
is the distortion of G. 

For this argument to work we need to take u such that J \Vu\\dz\ > 1 

for all 7 G r. In Q2 we take u(x,y) = (f^ f y (t)dt,y) and similarly for S7i . 
For each y the density f y (t) puts all the mass uniformly on the surviving 
intervals defining Qj and zero in the regions where we have no control on 
the distortion. We may say that for y G [2- Ni ,2~ Ni+N ], f y (t) = 2 Ni d~ i on 
the intervals defining ^2- Set Vit = 1 at all points where the distortion is 
bounded (the lower half of A and the top of A) . 

Now we have J \Vu\\dz\ > | J 7 Vu • dz\ > 1 so |Vu| is an admissible 
metric. 

Now we have that 

mod{G{t)) <C J J f*(x)D(x, y)dxdy + C 

<CJ2j 2 _ m / fy(x)Di(x)dxdy + C 
<C + C^ 2- m D i 2 2Ni d- 2i 2- m d i = C + C^ DidT 1 

i i 

So if Di increases slow enough the modulus will be bounded by a con- 
stant. 
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□ 

We will apply this theorem to an annulus where the distortion is con- 
trolled by D tk+1 < CAA tk+1 2 N e 5B ~ CA 2 A tk+1 for levels [2" nfc , 2- n *-i +1 ] 
to get Mod(G(A(z,r,2r))) > -gf- A — . The constant A k „ corresponds to the 

level r ~ 2~ nfe o . We are allowed to apply this theorem because the following 
relation holds 

oo 

r-^ n k Q +j- n k +i-l 

2^ A k)+i d n < oo (123) 

i=\ 

which is a lot stronger than the necessary relation 

oo 

r-^ "k +j- n k 

2^A ko+i d n <oo. (124) 



Inequality (|123[) holds because A k ~ Cne<T while — nfe-i ~ Cne<T • ^ n 



addition d ~ 2 n ~ 5 . 

We conclude this section with the following summary: if we run the stop- 
ping time algorithm for an annulus at scale 2~ nk o and obtain the surviving 
d-ary tree, the distortion of the modulus of image of the annulus is given by 
A 2 A k0 ^A^ e f. 

10 Putting it all together 

We are now in a position to give an outline of the proof of the main proba- 
bilistic estimate. Fix n and z. p n is fixed. We look of for p n , N n , b n and c n 



such that (11) holds 



If all of the following events hold for a large number of annuli 

• ti,k,t r i,k are small 

• There is an infinite d-ary surviving subtree in region Ri t \. 

• Li,k, $i,k,j are small, 

the sum of the moduli is big. Small sum of moduli implies a large number 
of failures for at least one of these rules. We have probabilistic estimates 
for each such failure. We will treat each rule separately and due to the 
decoupling we will be able to use the independence of variables corresponding 
to some of the scales. 
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We now proceed with the argument. 
Define the following indicator functions of events: 

i-1 

Xi (*) ■=Ylx(t i ,k<2 k - i In AO 

fc=0 

X2OO := X(*^i,i < l n Ai) 

fc>i 

00 min{iV n ,j-l} 

X4(»):= I] II x(^<2 fc -nn( J D„2 fc " J )) 

j=i+l fe=i+l 

XsW := X( L i,i < D n) 

X6(i) := II ^ k < 
k=i+l 

X7(i) '■= x(there is a good surviving d-ary subtree with A 2 Ai C0 = D n ) 

If nLi Xli?) = 1 then Mod(G{A(z, pnp^, 2~p n p i j2))) > cD~ 3 , where c is the 
constant from the modulus estimate and is independent of n and i. 

We divide all these events in three categories. Category I contains events 
of type 2, 5 and 7. Category II events are those of type 6 and all others are 
Category III. 

The event Mod(G(Ai)) < a n cD- 3 N n can happen for two reasons. 

Rl. either there is a set of 2a n N n of indices i for which all events of category 
II and III take place, but more than a n N n events of category I don't; 

R2. or no 2a n N n indices i exist for which all category II and III events 
take place. This can happen for two reasons. 

a. either for more than 1 — 4a n N n indices i category II events do 
not take place, 

b. or for more than 2a n N n indices i category III events do not take 
place. 

The following inequalities hold because of the independence of the XliP) 
for those particular values: 

. a n N n /3 



P(R1) < C2 N " ((1 + ln(A 1 ))e- ln2D "/ 2 ) Q " " (125) 
— ] +2 JV «(-^= 



+2 N -[^) + 2 N -(^=) (126) 



55 



The variable C n is the same which appears in lemma (17) and comes from 
the bound on the negative moments. 

We get a bound on P(R2.a) by following the argument in [AJKS09] using 
lemmata [17] and [T3j 

oJVnn/n n -l -(2+3o n /2) (l+26 n )jV„ 

whenever p n is chosen small enough to have 4p° n ^ 2 < 1/2 and D n is large 
enough to have C^-D" 1 p n ( 2+3a ™/ 2 ) < \ji (jj n \ s fa e con stant which appears 
in lemma (fl7|)). In addition a n and b n have to satisfy l+26 n = (l+a n /2)(l— 
4a n ). 

We present here the estimates for the sake of completion: 

P{R2.a)< £ E[l[j£(i)] (127) 

\S\=p ieS 

where p = 1 — a n N n and 5 is the set of indices where the events of type 6 
fail. We introduce the notation 

Xe(i,k) :=x{L i ,k<4 i - k D n ) 

The failure of a Xe(^) is due to some k, but it is possible that two different 
i's fail because of the same k. We have the bound: 

p r 

<Qx£(*)]<E E ^n^fe^i+^ ( 128 ) 

i£S r=l (i lr ..,Z r ) i=l 

where is the smallest i £ S such that X6(^i) = and level i\ + l\ causes 
X6(^i) to fail. More generally, is the smallest index i £ S larger than 
ij + lj and level ij+i + Zj+i causes X6(*i+i) to fail. 

The intervals + lj] cover the set S and this implies 

E ij > P - r 

3=1 

Applying lemma [17] we get: 

E\Xe&,ij + < CjiD- x p%- 1){ - 1+an) = C n L>-Vn (1+a " ) (4/^)^129) 
The terms on the right hand side of ( |128[ ) are independent so we get: 

^[II^«]<E(c , ^nVn (1+a " ) ) r E ( 4 Pn +a ") Eij (130) 
«G5 r=l (/i,...,Z r ) 



< p (l + a n /2) P ^^ njD - lp -( 2+ 3a n /2)j £ ( 4 ^/ 2 )^ (131) 
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where we have used the fact that Ylh — P ~ r - We now drop the constrains 
on lj to get: 



p / oo A r 

ElUx$(i)} < P ( ^ /2)P i: (CAVn-^f £ ( 4 ^" /2 ) (132) 

ieS r=l \2=1 / 

If we take p n small enough to have 4pn"^ 2 < 1/2 and D n large enough to 
have C n Z?-Vn (2+3an/2) < 1/2 we get: 

Mil X^)) < 2C„ J D- 1 p; (2+3a " /2 Vi 1+a " /2)p < ^ 1+a " /2)p (133) 

If we now take a n small enough such that 1 + 2b n = (l + a n /2)(l — 4a n ) 
we get 

E[U X 6 C «] < 2C nJ D- 1 p; (2+3a " /2 Vi 1+2b " )7V " (134) 

and 

P(i22.a) < 2 Af "2C„Z)-Vn (2+3an/2) ^ 1+26n)JVn (135) 

This concludes our bound on the failure of events of type 6. 

Failure of events of type 1 and 3 can be bounded as in [AJKS09J. We 
use a similar argument now to take care of events of type 4. The estimate 
we get for type 4 events also bounds the failure of type 1 and 3 events. 

The probability that more than a n N n events of type 4 fail is bounded 

by 

\S\=p ie5 

where p = a n N n and S is a subset of indices i of cardinality p. The letter C 
in the exponent of the indicator function refers to the complement/failure 
of the respective event. 

Each variable X4(0 = because of some bad event on level j. Although 
i's are different, the problematic levels j might not be; the same problematic 
level might cause problems for several i. This leads us to the following 
bound: 

p I oo \ / oo 

E\Y[xiii)]<Y, £ E xx(k,ji))-( E 

ieS r=l rn 1 +...+m r =p yj 1 =i 1 +mi J \j r =i r -\-m r 

min{N n ,j— 1} 

where X S, j) = \{ x(\k,j < HD^)) 

k=i 
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The bound should be read as follows: there can be between 1 and p problem- 
atic levels j. Say there are r problematic levels j. The total set of indices S 
will be partitioned in r subsets of cardinality mi, . . . , m r . If a level j causes 
problems for m indices i there must be an index i such that i + m = j. 
We now have the following: 

oo oo min{AT„J;-l} 

E x$W«)]< E E e -- a (^o^«)^^/4 

j;=ii+m; jl=il+mi k=i t 

If ii + mi < N n , this is bounded by 

E E+ E E 

jl=H+mi k=n ji=N n +l k=n 

Otherwise (if i\ + m/ > N n + 1) it is bounded by 

jl=il+mi k=i t 

Similarly, the second sum above is dominated by the term in which ji = 
N n + 1 and k = N n 

oo N n 

E E - Ce- ln2(D " 2JVn_li)p " g < Ce- ln2( - Dn2mi)p " q (137) 

ji=W„+l k=i t 

The first sum is dominated by 

N n Ji — 1 JV n 

E E^ E e- ln2 W _ii )^ (138) 

j[=i[+mi k=i t ji=k+mi 
< Ce~ ln2 ( Dn2mi ^ q < (7 e - ln2 ( D n)Pn 9 -21n21n(D n )m ; p- 9 -4 m 2 p -9 



Then we have 



E E (••■)■■■(■■■)< (140) 

r=l mi+...+m r =p 



P \ r 



^ e -rln 2 (D n )p n 9 e -2pln21n(D n )p n 9 | ^ Ce -Am 2 Pn q ^ ^-q 

(142) 



r=l \m=l 
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Here C is a universal constant and p n is small which implies the rightmost 
term is less than 2 r . We thus get: 

E[Y[X4 (*)] < Ce~ 2phl2ln{Dn)p " q - ln2{Dn)p - q (143) 
ieS 

So the probability that more than a n N n events of type 4 fail is less than: 
(j2 Nn e~ 2ctnNn ln2ln ( D ^)pn 9 ' — in 2 ] {D n )pn q (144) 
We conclude this section with the following inequality: 

P{Y^Mod{G{Ai)) < a n cD- 3 N n ) < C2 N " {{I + \n{D n ))e'^ D -' 2 ) " " (145) 



=i 



+2 N "(^) + 2 N -(- ( —) (146) 



+2 N "C n D- 1 p-( 2+3a "/Vp^ 2b ^ N fU7) 

+C2 Nn e ~ 2a " N ^ ln 2 ln(D„)p~ 9 -ln 2 (D n )p~ ^4g) 

as long as all the constants satisfy the necessary inequality given above. 
10.1 Final estimates 

We now proceed with an analysis of the order of magnitude of all the con- 
stants involved such that all the results hold and such that 

N n 

P(]T Mod(G(Ai)) < a n cD~ 3 N n ) < p n p£ +b ^ N " (149) 
i=i 

• a n is given to us by the size of the variance in the GFF. It's size is 
a n ~ n~ 37 . This forces us to take a n ~ a n /2 and b n ~ n~ 3 "' (earlier 
we had k in stead of n) . 

1/2 Z'yn^ 

• We also take C n < D n ' ■ Since C n ~ C^ eg we must have D n > 

^neg 

• The inequality Ap a n nl2 < 1/2. Since b n is less than a n it suffices to take 
4 P n" /2 < 1/2 which leads us to p n ~ 2- 10 ™ 37 . 

-1 -(2+3o„/2) 



We also need to have C„.D~W < 1/2. If we take D n > 

J neg 



CneT™ > this is automatically satisfied. 
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It suffices to have the following inequalities (these two terms dominate 
the others): 



(/-i \ a n N n /3 
im) - ~p^ +bn)Nn ^ 

2 N "C n D- 1 p^ 2+3a ^p ( n 1+2b ^ N " < PnP^ +bn)Nn (151) 
The second inequality can be rewritten as: 

2 N -C n D~ x p-^ a ^p h - N - < p n (152) 

For consistency we also need p n ~ 2- n ™ ~ pf 1 . . . p^ 1 ~ 2- 10 ( Ari +-+( ri - 1 ) 37Ar '- 

In order to obtain the desired modulus of continuity of /+ we need to 
have that a n N n D~ 3 is greater than a positive power of N n . 

If we set D n ~ C^™ 37 ™ the relations above become 

2 N nC -n s -t" 1 'a n N n /3 <• 2 -10(A r i+...+(n-l) 3 T JV n _i) 2 -10n 3 T {l+b n )N n 

C~^ Jn '' 2~ Nn < 2- 1 °( iV i+-+( n - 1 ) 37Ar ™-i) (154) 
given the way we chose p n . 

Since N n > a^D^ ~ n 37 C^ 7 " 7 it suffices to take: 

N n > lOiVi + . . . + 10(n - l) 37 iV„-i (155) 

This condition is satisfied by any sequence N n which is a power (larger 
than 2) ofC^ 37 " 7 . 

Take N n such that a n N n D- 3 ~ iV^ or iV n ~ Cnef ^ /e for 1/4 > 
e > (independent of n). Then 

Y,Mod(G(A t ))>Nt e . (156) 
i=i 

This implies the modulus of continuity is given by the relationship of 

W ( t ) ~ e - E^ 1 and t „ 2 " S?=i * 37 ^ . 

For t in the range [2" EIU , 2 ~ E^ 1 ^ N >] we get w(t) < w (2" ^ ^ ) < 
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We want this to be less than e ^ log « ^ . This implies we need to have 
the following: 



fci-ty < e a-SSLi* 37 ^ (157) 
which is equivalent to: 

n— 1 

(^o^ct^^^E^ 6 ^ ( 158 ) 



l-2e 



i=l 



n-1 



i=l i=l 

• The last terms on each side dominate the computation, so we need to 
have (n^Nn) 1 - 2 " < N^Zl- 

• If we set N n = e^ n ^ the last inequality becomes (1 — 2e)(/(n) + 
37 Inn) < (1 — e)f(n — 1). Rewriting this we get f(n) + 37 Inn < 
l + jbrJin-l). 

• If jTjy — > as n — > 00 (this is the case in our situation), we get the 
inequality above for any function f(x) for which < e for x large. 

• In our situation fix) ~ x'^ xl /e = e 37X ' 7 lnx / e and hence yj^j = 
^(7X 7_1 lnx + x 7_1 ). 

• If 7 < 1 for all x large enough y^j < e, which implies that we get the 



modulus of continuity e 



'(log I) 1 



Recall that we needed the following estimate in the proof of theorem [5] 
in section [3l 

Lemma 25. Let D(z) be the distortion inside D. Then 

E\ 



;[ / £><M] < 00 (160) 



Proof. In each Whitney square the distortion D is bounded by the constant 
C n (from the proof of lemma (l7[) for scales [2 _n ™+ 1 , 2 _n "]. Then we have 



£[ / DcL4] < / E[D]dA < C V2-""C n < C ^-""CjjfT" < 00 (161) 
Jo Jo „ „ 



□ 
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We conclude by the following summary: for the sequence t n = 2 — 4^, 7 < 
1, and for any e > 0, we consider the sequence n n = Yh=i * 37 ^m where 

Ni = Cneg . If we construct the random measure using this sequence we 
get a conformal welding and a modulus of continuity of e~( log 

For other sequences {t n } which converge to the critical value, we still get 
a conformal welding, but we can not prove the uniqueness by means of the 
removability theorems of Jones/Smirnov and Koskela/Nieminen. 

Finally, for t n — > 2, the best modulus of continuity for the welding map 
is going to be worse than Holder. 
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